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(57) ABSTRACT 

A navigation system comprises a cell-phone, a base cellular 
telephone site, and a webserver. Each is paired with a GPS 
receiver. The GPS receiver associated with the cell-phone is 
aided with information received from the cell-site and 
webserver that help reduce satellite search uncertainty. Time 
difference and/or frequency difference measurements are 
taken with data collected from what the cell-phone and its 
GPS receiver assume to be accurate time and frequency 
Correction information is used in post-processing of the 
velocity solutions computed by the cell-phone to arrive at 
more precise deteraiinations for the system. 

16 Claims, 6 Drawing Sheets 
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INFRASTRUCTURE-AIDING FOR 
SATELLITE NAVIGATION RECEIVER AND 
METHOD 

RELATED PATENT APPUCAnON 

U.S. patent application Ser. No. 09/687,044, filed Oct. 11, 
2000, and tided SATELLITE NAVIGATION RECEIVER 
AND METHOD is incorporated herein by reference. 

FIELD OF THE INVENTION 

The present invention relates to satellite-navigation 
receivers and systems, and more particularly to minimal- 
cost field devices that are assisted with information about 
precise time and frequency by a network to improve time- 
to-first-fix, accuracy, and manufacturing costs of the field 
devices. 

DESCRIPTION OF THE PRIOR ART 

Global positioning system (GPS) receivers use signals 
received from typically three or more earth-orbiting satel- 
lites to determine navigational data such as position and 
velocity. GPS signals are available worldwide at no cost and 
are now being routinely used to determine the location of 
automobiles to within one city block, or better. Dual- 
frequcncy carrier GPS receivers typically track a pair of 
radio caniers, LI and L2, associated with the GPS satellites 
to generate accumulated delta -range measurements (ADR) 
from P-code modulation on those carrier frequencies and at 
the same time track LI C/A-code to generate code phase 
measurements. Carrier frequency LI is allocated to 1575.42 
MHz and carrier frequency L2 is positioned at 1227.78 
MHz. Less expensive receivers tune only one carrier 
frequency, and therefore do not have adequate information 
to compute the local troposheric and ionospheric signal- 
propagation delays that appear as position errors. At such 
frequencies, radio carrier signals travel by line-of-sight. 
Thus buildings, mountains and the horizon can block 
reception, and multipath reflections can interfere with good 
reception. 

Each one of the constellation of GPS satellites in orbit 
about the earth transmits one of thirty-two unique identify- 
ing codes in a code-division multiple access (CDMA) 
arrangement. Such allows all of the many GPS satellites to 
transmit in spread spectrum mode at the same frequency, 
plus or minus a Doppler frequency shift of that frequency as 
results from the satellite's relative velocity. Particular sat- 
ellites are sorted out of a resulting jumble of signals and 
noise by correlating a 1023 "chip" code to one of the 
thirty-two pseudo random number (PRN) sequence codes 
that are pre assigned to individual GPS satellites. These 
codes are not necessarily being transmitted in phase with 
one another. Therefore, "finding" a GPS satellite initially 
involves searching various carrier frequencies, to account 
for Doppler frequency shift and local crystal oscillator 
inaccuracies. The searching also needs to find a code match, 
using 1023 different code phases and twenty or more pos- 
sible conelation code templates. 

The single largest uncertainty stems from the random 
frequencies possible from typical local oscillators at start- 
up. Therefore, the apparent-Doppler frequency is known 
only within wide search boundaries. Knowing the actual 
Doppler frequency is not much help, because the local 
oscillator can be so far off nominal on its own. 

From the user's standpoint, at least two operational char- 
acteristics of prior art GPS receivers interfere with complete 
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satisfaction. Such conventional receivers often quit working 
indoors because the buildings reduce the local signal field 
level to less than the receiver's maximum sensitivity. And, 
most receivers take a very long time to produce a position 

5 solution from a cold start. 

Intensive calculations in GPS receivers have necessitated 
high clock speeds and lots of expensive storage memory. 
These, in turn, demand expensive and comprehensive hard- 
ware. Manufacturers and users aUke would appreciate 

10 lighter and thinner navigation solutions that could use inex- 
pensive platforms or share pre-existing platforms for other 
applications. 

The Internet also represents a way for individual GPS 
receivers to monitor differential correction data in their areas 
and to off-load calculation-intensive tasks on regional web- 
servers that have high performance processors. 

SnapTrack, Inc., (San Jose, Calif.) is a commercial sup- 
plier of wireless assisted GPS (WAG) systems. Time, 

^ frequency, and approximate location data are extracted from 
a wireless network to assist GPS-signal processing in a 
navigation receiver. Such technology is described in a num- 
ber of United States Patents assigned to SnapTrack, includ- 
ing: U.S. Pat, Nos, 5,945,944; 5,663,734; 5,781,156; 5,825, 

^ 327; 5,831,574; 5,841,396; 5,812,087; 5,874,914; 5,884, 
214; etc. Also see, U.S. Pat. No. 6,078,290. 

SUMMARY OF THE INVENTION 

It is therefore an object of the present invention to provide 
30 a satellite -navigation receiver that can work indoors with 
extremely low signal-strength levels. 

It is another object of the present invention to provide a 
satellite-navigation receiver that produces position solutions 
rapidly after each cold start. 
35 It is a further object of the present invention to provide a 
satellite-navigation system that is inexpensive. 

It is a still fiirther object of the present invention to 
provide a satellite-navigation system that interfaces with the 
Internet. 

^ Briefly, a system embodiment of the present invention 
comprises a cell-phone, a base cellular telephone site, and a 
webserver. Each is paired with a GPS receiver. The GPS 
receiver associated with the cell-phone is aided with infor- 
mation received from the cell-site and webserver that help 
reduce satellite search uncertainty. Time difference and/or 
frequency difference measurements are taken with data 
collected from what the cell-phone and its GPS receiver 
assume to be accurate time and frequency. Correction infor- 
mation is used in post-processing of the velocity solutions 
computed by the cell-phone to arrive at more precise deter- 
minations for the system. 

An advantage of the present invention is that a system and 
method are provided that reduce the costs of user navigation 
equipment. 

Another advantage of the present invention is that a 
system and method are provided that improve sensitivity and 
time-to-first-fix enough for urban canyon and indoor use. 

A further advantage of the present invention is that a 
60 system and method are provided that makes higher level 
database services to be offered on the Internet which are 
related to real-time and historical user position fixes. 

These and other objects and advantages of the present 
invention will no doubt become obvious to those of ordinary 
65 skill in the art after having read the following detailed 
description of the preferred embodiments which are illus- 
trated in the various drawing figures. 
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IN THE DRAWINGS 

FIG. 1 is a functional block diagram of an infrastructure- 
assisted navigation system embodiment of the present inven- 
tion that includes a cellular telephone, a cellular-telephone 
base site, and a webserver; 

FIG. 2 is a diagram representing how the differences in 
two different clocks can be quantified as clock pulse accu- 
mulations in digital counters; 

FIG. 3 is a diagram representing a bias skew of a second 
event in the last one-thousandth part of a 1-PPS period; 

FIG. 4 is a diagram of the relationships between GPS 
events, cell-site events, and the various kinds of timing bias 
that are computed by embodiments of the present invention; 

FIG. 5 is a timing diagram representing the relationships 
between clocks and time periods in support of the math- 
ematical equations presented herein; 

FIG. 6 is a diagram representing a quantization of external 
cell-site clocks in a freqDiff observation; 

FIG. 7 is a diagram representing the physical geometry of 
a GPS satellite, cell-site, and cell-phone, and how timing/ 
frequency biases can affect position information; 

FIG. 8 is a timing diagram that represents each clock with 
its own system time and a corresponding clock bias from this 
time, the sum of these two produce a reference time used for 
making observations in a receiver, and the time -difference 
observations are taken in between these time-references; 

FIG. 9 is a chart of the frequency relationships amongst 
the cell-site, cell-phone, and GPS receiver 1-PPS outputs; 
and 

FIG. 10 is a functional block diagram of a second 
infrastructure -assisted navigation system embodiment of the 
present invention wherein no political or economic oppor- 
tunity existed to include a GPS receiver at a cellular- 
telephone base site, so a cell-phone and GPS receiver 
combination was permanently parked at a stationary location 
to serve instead. 

DETAILED DESCRIPTION OF THE 
PREFERRED EMBODIMENT 

Three hardware-based measurements can be used to syn- 
chronize the local clocks in two independent systems, e.g., 
(1) time-difference, (2) frequency difference, and (3) crystal 
temperature -versus-frequency models. Any or aU three can 
be used to synchronize clocks. In a timeDiff observable, if 
the time difference between time events from two different 
time sources is known, the time of an event at one time 
source can be used to compute the time of the event at the 
other time source. In a freqDiff observable, if the frequency 
difference between two clocks is known, then finding the 
frequency of one can be used to predict the frequency of the 
other. In a tempMeas observable, a temperature-frequency 
calibration model is used to predict the frequency of a local 
GPS clock by measuring the crystal temperature. If one of 
the time sources is a GPS receiver, other clocks can be 
related to GPS time. Such potentially allows devices all 
around the world to be time and frequency synchronized to 
a common and stable reference, for example, the GPS 
system atomic clocks. 

FIG. 1 illustrates a system 100, in an embodiment of the 
present invention, that includes a cell-site 102, at least one 
cell-phone 104, and a network server 106. The cell-site 102 
comprises a cellular transceiver 108, its reference clock 110, 
and a GPS receiver-112. The cell-phone 104 has a cellular- 
phone -service-band receiver 114 that is used to track the 
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frequency and phase of signals it receives from the cell-site 
102. The cell-phone 104 also includes its own GPS recciver- 
115. The network server 106 receives and distributes time 
information and other data from both the cell-site and 
5 cell-phone GPS receivers over a data network 116. A data- 
base 118 of positions and their time and frequency offsets 
from GPS time is collected for later reference. 

Passive network synchronization embodiments of the 
present invention use the cell -site GPS receiver to synchro- 
nizc the cell-site clocks to GPS time. Then both the cell- 
phone GPS receiver clock and the cell-phone clock can be 
synchronized to GPS time by the network server over a 
network infrastructure, e.g., the Internet. A typical GPS 
receiver generates a one-pulse-per-second (1-PPS) time ref- 
erence that can be compared and related with timeDiff 
circuitry to a 1-PPS clock derived from each of the cell -site 
and ceU-phone cellular-telephony radio carriers. Other 
observation intervals besides one second would also work 
the same way, but the mathematics are simpler when a one 
second time interval is used. 

20 

Embodiments of the present invention can use three 
different communication links to transfer time and frequency 
data. One links the GPS reference receivers 112 at the 
cell-site 102 with the network server 106. A second carries 

2^ software-API or packet-communication between the cell- 
phone 114 and GPS receiver 115 in the cell-phone, either 
directly in the cell-phone or indirectly through the network 
server. Sharing data directly between the receivers helps 
improve system performance. A third connection links the 
network server 106 and cell-phone to both the cell-phone 
and GPS receiver. The network server 106 can also be 
located at the cell-site 102 and communicate via the cellular 
telephony system transceivers 108 and 114. 
Time-difference -of- arrival measurements (TDOA), e.g., 

35 from a cell-site to a cell-phone, can be nsed for positioning 
inside a combined cell-phone and GPS receiver to improve 
positioning availability and accuracy. Time-stamps sent 
from the cell-site to the cell-phone can be combined with a 
priori knowledge of the cell-site position, the cell-site clock 

4Q bias and the cell-site drift from GPS time. Such time-stamps 
can be used to transfer time infonnation from the cell -site to 
the cell-phone with an accuracy affected only by the radio- 
signal propagation distance between the cell-phone and the 
cell-site. 

45 The corrective frequency offset of the cell-phone carrier- 
. frequency synthesizer loop can be used to estimate the GPS 
receiver frequency offeet. Such estimate is used to reduce the 
frequency uncertainty and thus reduce the timc-to -first-fix, 
improve receiver sensitivity, and/or reduce the size and 

50 power-consumption of the hardware. Digital numeric- 
controlled oscillator (NCO) and analog voltage-controlled 
oscillator (VCO) are both commonly used in such carrier- 
frequency synthesizer loops. 
The tempMeas measurement can be used to fiirther 

55 improve the frequency-transfer accuracy passing from a 
cell-phone to a GPS receiver. The freqDiff observation can 
be used alone in frequency transfers without any timeDiff 
observation. The timeDiff observation can also be used 
without the freqDiff observation for combined time and 

60 frequency transfers. The time transfer accuracy is preferably 
improved by reducing the timeDiff hardware quantization 
noise with software-based filters. Cell-phone observations of 
the cell-site frequency and phase can be used in dual-mode 
positioning system alternative embodiments of the present 

65 invention. 

The term "passive synchronization" as used herein means 
the time and frequency offsets of the clocks in the cell-site 
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and cell-phone are only observed with respect to a GPS 
receiver clock. Calibrations are made later by sending the 
naeasured offsets through the communication layers, e.g., 
using transfer control protocol/Internet protocol (TCP/IP) 
packets on the Internet. The real clocks in the cell-phone and 5 
cell-site are not controlled. The signals from these clocks are 
merely observed, and the corrective data obtained is not in 
any way fed back to the source to change the clocks. Rather, 
the data is used mathematically at the network server and in 
the GPS receiver software to assist in the GPS-sateUite 
search and then in the position calculation. 

FIG. 2 diagrams the time relationships between a GPS 
1-PPS per 202, a cell 1-PPS period 204, and an output 
counter value 206. A timeDiff observation can be made with 
a counter circuit that uses a GPS master clock for its input 
frequency. The number of internal clocks are counted 
between the two event sources, e.g., the internal 1-PPS 202 
and the external (or cell) 1-PPS 204. The result is periodi- 
cally output and the counter is restarted, such that for every 
second of lime there will be two coxmter outputs. 

Both output count values can be combined to estimate a 
key parameter, the bias between the two clock sources. Each 
restart at the external 1-PPS 204 produces a bias result. If the 
counter is instead restarted on the internal 1-PPS 202, the 
result is one minus the bias. The external 1-PPS (from the ^ 
cell-site or cell-phone) can be synchronized to the internal 
1-PPS (the GPS receiver) in one of at least two ways, e.g., 
with a smooth steered or an unsteered internal 1-PPS, 

Referring now to FIG. 3, it is assumed that a GPS 1-PPS 
timing edge is accurately slaved to GPS system time. In this 3Q 
way, the GPS 1-PPS timing edge represents "true" GPS 
time. A 1-PPS edge can be generated by a countdown circuit 
that is started by a 999th-millisecond interrupt derived from 
the GPS-integer second. An interrupt range of ±1 millisec- 
ond is preferred so that the 1-PPS can be adjusted to occur 35 
anywhere between msec-999 and msec-1 in the GPS-integer 
second. A msec-bias error is computed in the GPS receiver 
as the GPS -millisecond timer to "true" GPS-time error. If the 
GPS-receiver oscillator drifts from GPS time, such msec- 
bias error will change. The countdown value will also ^ 
change to reflect the drift. When the bias integrates more 
than ±0.5 msec, the GPS milhsecond counter is preferably 
adjusted to keep the circuit in range. 

The smooth-steered method embodiment of the present 
invention is used in GPS receivers at the cell-sites since the 45 
msec-bias error is known. Time is known in the GPS 
receiver and can be passively transferred to the cell-site by 
measuring any offiset between GPS time and a cell-site time 
reference. 

The time associated with the cell-site 1-PPS can be 50 
computed as, TimeCS equals timeGPS+biasCSG. Wherein, 
the variable biasCSG represents the time offset between the 
cell-site and GPS time. It is estimated using a function that 
uses either a snapshot or a filtered estimate using the 
time-series of the measurements Nil, Ni2, for all i. 55 

Hie unsteered internal 1-PPS method embodiment of the 
present invention does not use the additional countdown 
circuit. The internal 1-PPS is generated directly from the 
millisecond interrupt corresponding to the integer GPS 
second event. Such still requires that the millisecond inter- 60 
nipt that corresponds to the GPS integer second be periodi- 
cally shifted when the msec-bias error grows beyond a 
predetermined threshold. The time between intemal PPS 
events can be greater or less than one second. The step shifts 
that occur are brief, and do not adversely affect the receiver 65 
It knows when the event will occur, and the 1-msec shift can 
be accounted for in the computed residual msec-bias error. 
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The unsteered reference can be used by a GPS receiver in 
a cell-phone whenever the position is being requested from 
the cell-phone user interface. Such GPS receiver is typically 
not tracking satellites at the moment, and so the GPS 
msec-bias error will be unknown. The inverse problem 
exists in the cell-phone. GPS time is not known, and time is 
being transferred from the cell network to the ceU-phone, 

In operation, the time and frequency offset of a particular 
cell-site is obtained from a network server. The position of 
the cell-site is also useful in the solution calculations. The 
ofket between the cell-phone clock and the GPS receiver 
clock is measured. An estimate of the GPS receiver time and 
frequency offset is computed from GPS time, and is used to 
reduce the GPS satellite search range. 

The timeDiff observation is used in both the cell-phone 
and the cell-site to reduce the total number of unknowns in 
the system. The timeDiff measurement is defined herein as 
the difference in time between two time events generated by 
two different sources. Both clocks generate a signal with the 
same expected period based on the expected nominal fre- 
quency of each clock. The measurement is the number of 
clocks between consecutive events. 

Because all clocks have some firequency offset from the 
desired nominal frequency, the Bias is also constantly 
changing due to the effect of the frequency offset from 
nominal of both clocks. 

Suppose that "Freq** equals the actual physical frequency 
that is observed, it also equals baseFreq+drift (eq2). The 
"baseFreq" equals a desired frequency, e.g., the clock is 
designed to operate at this frequency but appears as "Freq" 
due to imperfections, aging, and temperature effects, the 
actual frequency appears. "Drift" equals frequency offset 
between the current frequency and the desired nominal 
frequency. 

Assume now that both clocks generate asynchronous 
one -second time event signal pulses with a specific time 
interval between. A one-second event can be generated with 
a counter that counts a specific number of reference clocks 
and then outputs a pulse when the desired count is reached. 
The one-second count number is the nominal crystal fre- 
quency. 

The actual time produced is equal to the number of clocks 
times the "true" period of each clock. Assuming a constant 
drift, the "true" period of each clock is the "ClockPeriod" 
variable, which equals the inverse of the frequency in Hertz. 

The actual time period produced, DT, equals baseFreq/ 
frcq. Inserting equation-2, 

DT = freq - driftlfreq (equation-3) 
= 1 - drift! freq 
= 1 ■»- deUaBias, 

where, deltaBias equals the change in the bias that occurs 
between the two events. 

The effect of equation-3 is that with each one-second 
event, accumulate an additional increment to the total bias. 

For the 1-PPS made with a steering clock when GPS time 
is known, the offeet can be estimated from GPS time and 
steer the clock to always be within one clock of the "true" 
time that was desired at the event. 

Alternatively, the 1-PPS can be left unsteered. The total 
bias is maintained in software and accounted for without a 
loss of synchronization when all calibration is in software. 

Extending to two clocks, bascGRfreq equals the expected 
the frequency of clock used in the GPS receiver (Hz). 
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"BaseCSfreq" equals the expected frequency of the clock 
used in the cell-site (Hz). "FrcqGR" equals the actual 
observed frequency of the GPS receiver (Hz). "FreqCS" 
equals the actual observed frequency of the cell-site (Hz). 
"BiasGRi" equals the bias of a GPS receiver's 1-PPS from 
"true'' GPS time at time Ti. "BiasCSi" equals the bias of a 
cell-site's 1-PPS from "true" GPS time at time Ti. "Delta- 
BiasGR12" equals the change in bias in the GPS receiver 
from tl to t2. "DeltaBiasCS12" equals the change in bias in 
the cell-site from tl to t2. 

A reference clock and an unknown clock are defined, then 
the GPS clock is chosen to be the reference. It can synchro- 
nized to GPS by tracking the satellites. 

The cell-site is defined to be the unknown clock that is to 
be synchronized to GPS time. 

The "true" time between the GPS and cell-site events is 
defined. TDij equals "true" time between two consecutive 
events in the i-th time where the source of the newest event 
is the j-th source. J=l is the cell-site and j=2 equals the GPS 
receiver. 

FIG. 4 represents the relationships between the "true" 
GPS 1-PPS events, the actual locations of the 1-PPS events, 
and the "true" tDiff observations. These observations can be 
made. 



73,030 Bl 

8 

FIG. 5 represents a redrawing of the timeline and allows 
the effects of the quantization enor to be seen. Nil equals 
three because the first cell-site PPS is reported on MCLK3 
due to the clock synchronization. The "true" measurement 
5 TDll is equal to the measured count (TDllm) minus the 
quantization Ql. 

The "tme" time difference is the "true" number of clocks 
(integer plus fraction, which is the quantization) times the 
period of the clock. Assume the frequency was constant over 
10 the interval, then the period is the inverse of the "true" 
frequency. The "truth" model for the ideal measurement is 
shown below. 

Hie "true" time differences are defined as, 

15 TDll equals {Nll-Ql)JfreqGR 

TDM equals (m2+Ql)/freqGR 
TDll equals (N2l-Q2)/freqGR, 

20 The measured counts are integers. The exact frequency of 
the clock used to count is only an estimate. The measure- 
ment used is the observed number of counts, times the 
predicted period which is the inverse of the expected fre- 
quency. 

25 The measurements are defined by. 



BiasGR2 - biasGRl + deltaBiasGRlZ 

- biasORl - driftOR/ftcqOR 
BiasCS2 - biasCSl + deltaBiasCS12 

- biasCSl - diiftCS/freqCS, 
TDll = (1 + biasCSl) - (1 + biasORl) 

- biasCSl - biasGRl 

TD12 - (1 + biasGRl + (1 - driftGR/ficqOR)) - (1 + biasCSl) 

- 1 + biasGRl - driftGR/frcqGR - biasCSl 

- 1 + biasGRl + deltaBia5GR12 - biasCSl 

- 1 + biasGR2 - biasCSl 

- 1 - (biasCSl - biasGR2) 

TD21 - (1 + biasCSl + (1 - driftCS/fircqCS)) - (1 + biasGRl + (1- 
driflGR/freqGR)) 

- (2 + biasCSl + deltaBiasCS12) - (2 + biasGRl + delta BiasGRl 2) 
» bia5CS2 - biasQR2 

1D22 - (1 + biasGRl + (1 - driftGR/freqGR) + (1 - driftGR/frcqGR)) - 
(1 + biasCSl + (1 - driftCS/frcqCS)) 

= 1 + biasGRl - driftOR/freqGR - driftGR/frcqGR - (biasCSl - 
driflCS/frcqCS) 

» 1 + biasORl + deltaBia50R12 + dcltaBiasGR23 - 
(biasCSl + deItaBiasCS12) 

- 1 - (biasCS2 - biasGR3). 



Summarizing, if cell -site is the source at time t«i, then j=l 
and, 

TDil'bUtsCSit^O-biasGRit^HX (cq4) 

if GR is the source, at time t=I, then j=2 and 

TDa-\-({biasCS{t-i)-biasGR{t^i+l)), (eqS) 

The quantization error is harder to remove when only the 
timeDiff observable is available, and this motivates gather- 
ing the freqDiff measurement. Such allows the quantization 
noise to be reduced by filtering, and higher frequency clocks 
are not needed. Avoiding higher frequency clocks allows a 
system operating power savings to be realized, and is very 
important in battery-operated portable devices. 

Two important parameters can be estimated from the 
timeDiff observables, the bias from GPS time of the cell-site, 
and the frequency drift of the cell-site. The frequency drift 
estimate is then correlated with the bias estimate, and no 
new informatioD is needed to filter the cell-site bias esti- 
mates. 



TDUm^mXIfnqGRhat'^lfniqGRhat 
TD\.lm«m.llfrtqGRhaU'VfreqGRhat 
30 TD2lm=N2 llfre{(GRhatm2fFrtqGRhat, 

where, 

FreqGRhatoestimated GPS receiver frequency 
FreqGRhat^baseGRfreq+driftGRhat 
DriftGRhat=estimated frequency offset from nominal 

made by the GPS receiver. 
The "tme" time difference model defines, 

TDlUbiasCSUbiasGRl 

40 

TDM'^l-UasCSUbiasGRl 

TDZUbiasCSZ-biasGKl. 

45 The numeric index on the bias parameters represents their 
time-tags. The first numeric indicator on the TD is the time, 
the second is the source used to trigger the observation. 

The first measurement to find the error model for the 
estimate of biasCSl starts by defining the error to be the 
50 "true" one minus the estimated parameter. 

The estimation error is defined as "errorBiasCSll", and is 
the estimation error of the cell-site bias in the l"" second 
using the cell-site PPS as the trigger. 



ErrorBiasCSll 

» biasCSl = biasCSlhat 

= CTDll + biasGRl) - CTDllm + biasGRlhat) 

- ((Nil - Ql)/frcqGR + biasGRl) - 
(Nll/freqGRhat + biasGRlhat) 

The auxiliary variables are defined, 
EnorBiasGR - biasOR - biasGRhat 
EnorDriftGR » freqGR - freqGRhat 

- (baseGR&eq + driflGR) - (baseGRfreq + driftGRhat) 

- driftOR - driflGRhat. 
Continuing, 

crrorBiasCSll 

- Nil * (1/freqQRl - l/frcqORlhat) - 
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-continiied 



Ql/&cqGRl + crrorBiasGRl 

- Nil " (freqGRlhat - ercqGRl)/(£rcqGRl * frcqGRlhat) - 
Ql/&eqGRl + crrorBiasGRl 

- -Nil * crrorDriftGRl/(fircqGRl * freqORhatl) - 
Ql/freqGRl + crrorBiasGRl . 

Assuming Ql equals 1, crrorDriftGR equals .01, driflGR - 0, 
and bascGRfreq = 27.456 MHz, then the error is, 

errorBiasCSll - TimeDiOVue • (36.42191142 ns - 
36.42191142 ns) - 1 '36.42191142 ns 

a 36 OS 

- one clock of the master clocL 



This is the worst-case accuracy error, which is about 10.91 
meters. 

The first term can be ignored because the GPS drift error 
is much less thao the total frequency. In general, the GPS 
receiver can measure frequency to less than 1 Hz at LI. The 
error at the nominal frequency is 1/(1575.42/27.456), which 
equals 1/57.4 and 0.0174 Hz at bascGRfreq. Since Nil can 
be close to freqOR, this breaks down to 0.0174/freqGR 
which is 0.63 nanosecond equals 0.19 meter. Such term can 
be ignored when the drift error is small. 

The second term is the larger error source. Worst case, Ql 
is close to one clock, and Ql/baseGRfreq equals 36 nano- 
seconds or 10.91 meters. If both edges of the master clock 
are used, an accuracy of eighteen nanoseconds can be 
realized. Both such edges of the reference clock are used to 
drive the counter between the consecutive 1-PPS events. The 
resolution is improved and the worst-case quantization is 
_/baseGRfreq equals 18 ns equals 5.46 meters. 

The above example gives the theoretical worst-case accu- 
racy for TDiff based on only the frequency of the reference 
clock, which is 27.456 MHz. 

A second measurement in the first interval provides a 
second estimate of the cell-site bias, but with a different GR 
bias, 



ErrorBiasCS12 

- biasCSn - biasCSnhat 

« (1 - TD12 ■»- biasGR2) - (1 - TD12m + biasGR2hat) 

- (-(N12 ^- Ql)/freqGR2 + bia6GR2) - (- 
N12/frcqGR2hat + bia5GR2hat) 

= -N12 (l/freqGR2 - l/freqGR2hat) - 
Ql/freqGR2 + errorBiflsGR2 

= -N12 CfreqGR2hal - freqGR2)/(£reqGR2 * freqGR2hat) - 

Ql/freqGR2 + errorBiasGRl 

« N12 (crrorDriftGR2)/(eTeqGR2 * frcqGR2hat) - 
01/freqGR2 + errorBiasGR2. 



The GPS drift estimation error is small so the first term can 
be ignored. Two estimates are obtained for the cell-site bias 
in one second with different noise errors. 
The estimated measurement models are, 



ErrorBiasCSVl^QVfrtqGRl^rrorBiasGia^ 
m2(errorDriftGR2)/(fnqGJa*fKqGRZhat) 

ErrorBlasCSli'-QllfreqGRl-^rrorBiasGia- 
Kl\{errorDnftGia)j(freqGRl''fnqGP2haty 

5 

Using both edges of the 27.456 MHz GPS clock, the 
quantization error is a uniform random variable with a 
maximum error of l/(2*27.456e6), and equals 18 ns or 5,46 
m. 

10 The GPS errors can also be modeled, 

et7orBiasGR-N{0yiigm oBias 1) 
sigmaBias^TDOP *sigmaPR 

sigmaPH = (l-2 meters ) / speed -of- Eght 
= 3-6 nanoseconds. 



20 errorDriftGR~N{Q^igmaDnft1) 

stgmaDrifU'TDOP*sigmaPRR 

stgmaPRR=OM m/s=0.003 nanoseconds/second. 

25 where, N(mean,variance) is a Gaussian distribution. TDOP 
equals time-dilution -of -precision from the least-squares 
solution matrix. The numbers given above for the standard 
deviations for the pseudo-range (PR) and pseudorange-rate 
(PRR) are for conventional tracking with code and carrier 
30 tracking loops. 

An estimate of the cell-site clock drift can be formed 
using timeDiff data. Two time^difference measurements can 
be added together to span two external PPS events from a 
cell-site. The sum of these is an estimate of the time of the 
35 cell-site one-second event, 

TDJ2 + TD2J = (i- biases i + biasGRl) -t- {biasCS2 - biasCRl) = 
1 + biasCS2 - biasCSI = 1 -t- biasCSJ ^ 
^ deiioBiasCSn - biasCSJ = I - dripCSIlf freqCS, 

where, DriftCS12 equals drift in cell-site clock in the 
interval between the cell-site PPS events. 
45 The frequency can be assumed to be constant over the 
interval so that freqCSl equals freqCS2, 

Substituting, TDi2 + TD2! = 
50 I - driftCSI2/ freqCS = I - drifiCSJ2/ {baseCSfreq + driftCSI2). 



Simplifying the derivation, TD=l-d/(fwjc-Ki) 
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55 TD*{6a5C+0-(basc+d)-(/ 

BiasCSllhat^VfreqGRViat^UasGRVtat TD^base^^base 

BiasCS21hatml-N2lffreqGR2hat't-biasGR2hat •d^base{1-TD) 

BiasCSl2hat-N12/fntqGR2haUbiasGR2haL d~base{l/TD-l). 

Where Nil, N21, and N12 are the counts from the hardware Expanding back to the usual notation, DriftCS12 equals 

measurement and freqGRhat and biasGRhat are calculated baseCSfreq* (1/(TD12-TD21)-1). 

using the GPS receiver. Substituting the definitions of the timeDiff measurements 

The error models of the cell-site clock biases are, "^es the equation for the "true" driftCS, driftCS12 equals 

65 baseCSfreq*(freqGR2/(N12+N21+Ql-Q2)-l). To define 

ErmrBiasCsn^Ql/freqGRi+emrBiasGRi-Nn(errorDriftGMi)/ the driftCS12 estimate, driftCS12hat equals baseCSfreq* 

{fr^qGR\*ff^qGRihaf) (freqGR2hat/(N12+N21)-l). 
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In order to simplify the error in the drift estimate, return 
to the format of the third equation in the simplified notation 
above. Also multiply both sides by frcqGR 

TD*{base-i-d)~base 

The "true" model, 

(CAn2+ei)+(/V21-e2)) • (baseCSfmq^driftCSi)^ 
baseCSfreq*freqGRl. 

Applying the same to the measurement model (except that 
multiply both sides by freqGR2hat), 

{im+ha I) • {baseCSfrsqHiriftCSVuU)''baseCSfKq*freqGF2hat. 

Subtracting the measurement model from the "truth' 
model, 

{Ql-Q2y{baseCSfreq+driftCSl)+(N2UN2iyidriftCSl- 
dnfiCSlhat)'baseCSfreq(frvqGRl-freqGmhat) 

Applying the following definitions, 

freqGR-baseGRfreq^-driftGR 
freqRhat'^bascGRfreq-^driftGRhat 
ErrorDrifiCS'driftCS-drifiCShat 
ErrorDrifiGRmdriftGR-driftGRhat, 

{Ql~Q2y{baseCSfrvq+drifiCSl)^mi+N2iy(drifK:S~driftCS- 
hat)''baseCSfnq*(baseGRfreq+driftGR2-baseGRfivq- 
dnfiGRlhat) 

{Ql-Qiy{baseCSfmq+driftCSl)+{mUN2iy{ertvrDriftCSl)' 
baseCSfnq*envrDriftGR2 

Rearranging to isolate the errorDriftCS term, 

(N2UNliy€rmrDriftCSX'baseCSfTvqX*errvrDriftGRl-igi- 
QZyfKqCSi. 

The previous oell-site bias estimate is protected forward ^ 
in time between the two cell-site PPS events. 

A new parameter is defined that is the projection of the 
previous biasCS estimate with the drift estimate to generate 
a reference trajectory to filter the new biasCS estimate 
against. 

The "truth" model. 
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BiasCS2Mmus~biasCSn-driftCSl/{baseCSfreq4'dnfiCSl). 

The measurement model, 50 

BiasCS2MinusHat--'biasCSi2hat-driftCSlhat/(baseCSfreq+ 
dnftCSlhat). 

Multiplying both sides by their respective estimated cell- 55 
site total frequency, 

Bia3CSminus*{baseCSfrt(iiiinftCSl)''biasCSll*(paseCSfTzq¥ 
drifiCSVj-dnftCSl 



BiasCSlMinusHat *(baseCSfKg+driftCSlhat)~ 

bi(uCSl2hat*(baseCSfreqi-driftCSlhat)-driftCSlhat 

Forming the "true'' minus the measured models, 

(BiasCSlMinus-BiasCSlMinusHatybaseCSfreq-^ 

SiasCS2Mitms*driftCSl-BiasCS2AfumsHat *driftCSlhatm(bi- 
asCSll-biasCSnhatybaseCSfreqi-biasCSn *dnfiCSl- 
btasCSl2hat*dHftCSlhat-(drifiCSl-driftCSUai). 



60 



The projection error is defined as, 

ErrorBiasCSZminus equals biasCS2minus-Bi<uCS2MinusHai 

Using the convention that hat equals truth-error, 

envrBiasCS2minus*b(ucCSfreq^bicisCS2Minuj*driftCSl- 
(biasCS2Minus-errorBiasCS2rmrmsy(driftCSl- 
erwrDriftCSl)" 
. erwrBiasCSU *baseCSfKq+biasCS12 " driftCSl- 
{biasCSn~crrorBiasCSny{dnftCSi-erwrDnftCS\)- 
ermrDriftCSl, 

envrBiasCS2minus*baseCSfreq+€rrvrBiasCS2minus*dnftCSi+ 
{biasCS2\finus-errorBiasCS27mrms) "errorDrifiCSi" 
envrBiasCSn*baseCSfKq+erFvrBiasCSi2 *driftCSl ■»- 
{biasCS12-errvrBiasCSl2) *errorDriftCSl- 
crwrDriftCSi. 

For the convention that error equals "lruth**-hat, 
then tmtb-erior equals hat errorBiasCSl 

mh\us*{baseCSfreq+drifU:Sl)-i-biasCS2MitwsHat*errorDriftCSl^ 

errorBiasCS12*{baseCSJreq+driftCSl)+ 
biasSCnhat*errorDriftCSl-e7TorDrifiCS\ 

errorBiasCS2minus*frtqCSl~ 

envrBiasCS12*freqCSl+(biasSC12hat~ 

biasCS2MinusHatyerrorDrifiCSl-errorDriftCSl 

Dividing both sides by freqCS and rearranging, 

ErrorBiasCS2minusm 

(iaotBmsCS12+(biasSCi2hat-biasCS2\fmiuHat-iyer- 
rorDriftCSVfreqCSl. 

Plugging in the definition, 

BiasCSlMinusHat equals biasCSHhat- 
driflCSihatl{baseCSfreq+dnfiCSihat). 
And rearranging, 

biasCSHhat-BiasCSlhatMinm equals 
driftCSViatKbaseCSfreqi-dnftCSihat) 
ErrorBiasCS2minus<' 

errorBiasCS 12+(- 1 +drifiCSlhat/{baseCSfreq+ 
driftCSlhat)yerrorDrifiCSllfreqCSl, 

Simplifying the term in the brackets, 

Using more simple notation: -l+d/(b+d) equals (-b- 

d-KJ)/(b+d)=-b/(b+d) 

errorBiasCSlminus^ 

erro rBiasCS 12- (fcaj<^Ji/«g/(baseC^//«g+ 
dHfiCS\hat)yerrorDrifiCSVfreqCS\. 
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Substituting that: baseCSfrcq/(baseCSfreq-KiriftCSlhat) 

equals (N12+N21)/freqGR21iat, 
errorBiasCSlmirmS'' 

moiBiasCSl2-{(^+J^l)/fn:qGR2hat)*ermrDnftCSiJfnqCSi. 
Plugging in the definition of 
errorBiasCSll' 

-QlJfreqGR2+errorBiasGm+m.2(erTorDrifiGRiy 
(fivqR2*freqGRZhat). 



Using the result in deriving the error in the drift estimate 
after dividing both sides by freqGR2hat, 15 



{(N2UN21)/fKqGR2hat)*envrDnftCS' 

'{Ql'QZYfniqCS\JfrcqGRlhat-^baseCSfreq*crrorDnftGRlf 
fr^qGRV\au 

Inserting these two equations into the working equation, 

envrBiasCS2mimts°-Ql/freqGR2+erTvrBiasGR2+ 
mi{erwrDnfiGR 

2)f(fKqGR2*ffvqGR2hat)+(gi- 

Q2)lfreqGKihax-baseCSfrtq *errvrDriftGR2}(freqCS 

l*fntqGR2kat), errorBiasCS2mimts~Ql *(llfreqGRlhat~ 

VfreqGR2)-Q2/freqGRhat+envrBiasGR2+ 

{errorDriftGRl/fnqGRlhat) * (NU/fmqGRZ- 

baseCSfreq/freqCSl), 

errorBiasCSlmimts~Ql *(freqGR2~freqGR2hat)/ 
{freqGR2hat*freqGR 

2)-Q2/freqGRhat+ermrBiasGRZ-t- 

{erTorDHfiGR2/freqGR2hat) * (NU/freqGR2- 

baseCSfreq/freqCSl), enorBiasCSlmiimsm 

Ql *errorFrtqGR2lifTtqGR2hat*frtqGPl)- 

QVfreqGRhat-^errorBiasGRl+ierrorDrifiGRl/freqGRl 

hat) * Qm/firqGRZ'bojeCSfrcq/fn^qCSl). 

Remembering the definitioD that, 

baseCSfreqIfrtqCSl equals {Nl2^N2l)lfreqGR2^{Q\' 
Q2)/freqGR2. 

Simplifying the last term in parenthesis as, 

{m2/freqGR2-baseCSfKq/fnqCSl)- 
Nl2/freqGR2-(liNl2+N2i)/frsqGR2'¥(Qi-Q2yfTvqGR2)- 
-N2VfreqGR2-(Ql-Q2)/fKqGm 
erre>rBiasCS2minus- 

Ql*errorDriftGR2((fKqGR2hat*fnqGR2)- 
Q2/freqGRhat-i-errorBiasGR2+(ermrDriflGR2/freqCR2 
hat) • {~N2lffr«qGR2-(S3l-Q2yfreqGR2). 



errorBiasCS2mmus» 

QVcrrvrDriftGR2/(fivqGR2hat*fif;qGR2)- 

Q2/freqGRhat+ 

cnorBiasORZ- 

N2VermrDrtftGR2/(fnqGR2hat*frt:qGR2)- 

(Ql-Q2yem)rDr{fiGR2f(fnqGRZhat*fnqGR2) 

errorBiasCSlminusm 

Ql'errorDriftGRV(fKqGR2hat*fKqGR2y 

Q2/freqGRhat+ 

enorBiasGR2- 

N2VenvrDriftGR2/(fKqGR2hat*frr:qGR2y 

{QUQ2)*crrorDnftGR2i(frcqGRZhat*frtqGRl) 

errorBiasCSlmmus' 

-Q2/freqGR2hat-t- 

crrorBiasGR2- 

N2V€rrorDrifiGR2/(fnqGR2hat*freqGR 

2)Q2^eTrorDriftGR2lifrcqGR2hat*intqGR2) 

errorBiasCS2minus^ 

Q2{-fnqGR2'hermrDriftGR2)/(freqGR2hat*freqGR2)+ 
errorBiasGRl- 

N21 • erwrDnftGR2/(fKqGR2hat*f7eqGR2) 
errorBiasCS2tninuS' 

~Q2{fr^qGR2~errorDriftGRl)l{freqGRlhat*freqGFl)+ 
erwrBiasGRl- 

N21*envrDnftGR2/(fnqGR2hat*freqGR2) 
errorBiasCSZminus- 

-Q2(frKqGR2hai)!(freqGR2hat*freqGR2)+ 
errorBiasGR2~ 

N2\' errvrDrifiGR2/(freqGRZhat*freqGR2) 
errorBiasCS2minus^ 
-Q2/freqGR2+errorBiasGR2-N 
2VerrvrDriftGR2/(fKqGR2hat*freqGR2) 



55 This equals the equation derived above for errorBi- 
asCS21. 

Quantization and GPS receiver errors will appear in the 
final filtered estimate of the cell-site, A complimentary filter 
60 is used to combine the two estimates of the same parameter 
into a single filtered estimate. The complimentary filter has 
the form, FilteredEstimate equals (1-a) 
*previousFilteredEstimate+a*ncwUnfilteredMeasurement, 
where a <=1.0. 

65 If both the previous filtered estimate and the new unfil- 
tercd measurement are estimates of the same parameter, then 
the filtered estimate is imbiased, 
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envrBiasCS2mmasmQl*(lffreqGR2hat'VfreqGm)- 



FilteredEsrimate equals {1 -a)*iX 4-Ni)4-a*{X ^ N2) = 
((l-a + a)*X)+(fl*/V2 + (l -a)*NJ)) = 

X ^i{i-a)*Nl + 0*^2). 5 

If the gain parameter "a" is small, then the noise on the 
previous filtered estimate is preferably not filtered. If there 
is an error in the propagation from a previous time to a lo 
current time, the error is placed directly into the filtered 
estimate. 

An optimal filter has a gain of 0.5 when the propagation 
error is on the same order of magnitude as the noise in a new 
raw measurement. In this way, some attenuation of the 15 
propagation error is obtained. 

Applying the compfimentary filter, BiasCS12filtered«(l- 
a)*biasCS12minusHat+a*biasCS12hat. The error then has 
the form, Error-(l-a)*(error in biasCS12minusHat)+a* 
(error in biasCS12hat). 20 

It is not possible to reduce any of the error sources 
because the predicted bias is equal to the measured quantity 
and the error in both terms is the same. 

An earher estimate of the cell-site bias can be used other 
than that allowed by the GPS observation. Variable error- 25 
BiasCSll is used instead of errorBiasCS12. 

Starting with the "truth" model, BiasCS2Minus- 
biasCSll-driftCSl/(baseCSfreq+driftCSl). 

The measurement model, BiasCS2MinusHat- 
biasCSllhat-driftCSlhat/(baseCSfreq+driflCSlhat). 30 

Arriving at, 

errorBiasCSlminus'' 

cnOTBUisCSn-(baseCSfreq}(baseCSfreq+driftCSlhat 
))*ermrDnftCSVfreqCSl. 

Substituting that, 

40 

bastCSfreq}{baseCSfrtq+driftCS\.hat) equals 
{NU-^Nl lyfreqGRZhat, 

erwrBiasCS2minus^rwrBiasCSl2-((Nl2+NZl)/freqGR2hat 

45 

) 'errorDrifiCS/freqCSl. 

Plugging in the result, 

erwrBiasCSlU- 

QVfre(fiRUcrrorBiasGR\-mX{ermrDriftGRl)l(frc(fiR ^° 
\.*frtqGRVxat). 

And using the result to derive the error in the drift 
estimate after dividing both sides by fireqGR2hat, 55 

-(Ql-Q2yfreqCSVfreqGR2hat+baseCSfreq*errorDriftGR2/ 
fn:qGR2hat. 

60 

Inserting these two equations into the working equation, 

envrBiasCSlmirms'^Ql/freqGRl+ermrBiasGRi-N 
n{crrorDriftGRl)}ifrtqGFX *fKqGRlhat)+(Ql- 
QiyfritqGRlhat-baseCSfreq *errorDnftGR2/(freqCS 
l*fTtqGF2hat), 



Qlffrtq^R2hat^rroTBiasGR\.-mi{erroTDriftGRl)lifrtqGR 

l*fKqGRViat)-baseCSfKq*errvrDnftGR2/(fKqCS*/KqGRZhat\ 

erwrBiasCSlminus~QX*{fre(fiR\-freqGRlhat)}(JregGRl 

h&l*(TcqGKi)-Q2/freqGR2hati-€rrorBiasGRi-N 

n{errorDriftGRi)/(freqGRi *freqGR\hat)- 

baseCSf req ''errorDriflGR2/(/>-e^Sl "freqGRIhat), 

errors iasCS2mmus=Ql''(freqGRl~ 

fKqGR2+ermrDrifiGRiy(fnqGR2hat*frvqGRl)- 

Q2/fnqGR2hat+ermrBiasGRl-mi (errorDriftGRl)f(frKqGR 

X^frtqGR\}mt)-baseCSfreq*errorDnftGRltifreqCSX*freqGRlhai). 

Thus, the difference of biasCS2hat-biasCS2minusHat, 

ErrorBiasCS2 i-errorBiasCSlminus- 

-'Q2/fr€qGR2+erFvrBiasGR2-iai{errorDriftGR2)/(fKqGR2* 

£rcqGR2/wr)-ei * (freqGRl- 

freqGR2+errorDn/tGR2)/(freqGR2hat'*freqGRl)+ 

Q2/freqGR2hat-em}rBiasGRUNn(errorDrtfiGRl)/(fnqGR 

l*freqGRlhat)+baseCSfreq*errorDnfiGR2f(freqCS 

X*frtqGRlhat)yQl*{\tfrrqGRlhai-VfreqGR2y 

Ql*{bas€GRfreq+driftGRi-bas€GRfreq- 

^mGB2+errorDriftGR2)/(fKqGR2hat*frriqGRl)+ 

errorBiasGR2-errorBiasGR I + 

mi{crrorDriftGRl){{freqGRVfrcqGR\hat)- 

N2 \{errorDriftGR2)lifreqGF2 *freqGR2hat)+ 

baseCSfreq *crrorDnftGR2{(frKqCSl *freqGR2bat)^ 

Q2 *errorDriftGR2}(freqGR2hat*freqGRl)' 

Ql*{driftGRl-drifiGR2+ermrDnftGR2)/(freqGR2 

\at*(rcqGKiyi-errorBiasGR2-errorBiasGRl+ 

mi{errorDnftGRl){(freqGR\*freqGR\hat)- 

N2 l(errvrDriftGR2y(fr^qGR2 *fTvqGR2hat}+ 

baseCSfreq*errorDnftCRll(frtqCSt ^fretfiKlhat). 

It appears that this may actually be smaller than the error 
derived earlier. However, this result is inconclusive and 
needs more investigation. 

Counting the number of clock pulses from the clock witb 
unknown frequency error in a known time interval generated 
by the reference clock, forms the freqDiff observable. The 
frequency is computed as the number of clocks divided by 
the "true" time interval. Even though there is enough infor- 
mation to make this estimation with the timeDiff measure- 
ment data, it turns out that the error in the estimate is 
perfectly correlated with the bias estimates. The same 
observations, with their associated quantization sources, are 
used in estimating the bias and drift. Thus, freqDiflf is used 
to form an estimate of drift that is statistically independent 
from the bias estimates. 
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Embodiments of the present invention count the clock 
pulses with unknown frequency error directly. In contrast, 
the timeDiff circuits count pulses from the reference clock. 

The quantization error is a uniform random variable 
between (0,1) and it is not zero mean. To reduce the 
parameter, the data must be pre -whitened to produce a 
zero-mean measurement error, and conventional techniques 
can be used in parameter estimating with a quantized noise 
souce. 

Once the noise souce becomes a zero-mean process, it 
may be reduceable with good filtering. However, in estimat- 
ing a parameter that has a non-constant slope, the previous 
filtered bias must be propagated to the current time. This is 
needed to assimilate the new measurement with the previous 
data in an xmbiased manner. The difference between the 
average bias and the new data is scaled and then added to the 
propagated estimate to finish the update. 

The 1^^ order low-pass filter, 

Bfiltered(K) equals PropagatedBtas(K)+gain 
{K)\biasMeas{K)~PropagatedBias{K)), 

where, 

PropagatedBias(K) equals Bfiltered(K-l)-drift(K) 
*deltaTime(K) 

gain(K) equals 1/m, where 
m equals k if k<Kmax 
m equals Kmax if k>Kmax 
deltaTime(K) equals time between t«K and t=K-l. 
Note, if Kmax equals infinity, then the filter output is the 
average of all the data. 

The propagation error from the timeDiff based-drift esti- 
mate has the same error as the new measurement. The 
propagation error is not filtered, and cannot use the timeDiff 
based-drift estimate in the projection phase. 

Embodiments of the present invention directly count the 
unknown frequency error clock pulses. The timeDiff circuit 
count pulses from the reference clock. 

The freqDiff measurement is the number of external 
clocks counted in an interval defined with an integer number 
of internal clocks. An estimate of frequency offset of the 
internal clock can measure the frequency of the external 
clock. 

The "true" measurement Is, ClockCount equals 
f reqCS * DeltaTimeTrue. 

The observed count will be in error by the effects of the 
quantization of the external clock against the internal master 
clock at the time the count is started and finished. 

The observation model can be derived from FIG. 6, in 
which the circuit counts four external clocks. However, the 
real number is a number that contains a fractional compo- 
nent equal to (Q1-Q2) where Ql and Q2 have units of 
counts and are continuous random variables between 0 and 
1 count. 

The "true** integer plus fractional number of counts, 
ClockCountTrue equals freqCS*DeltaTimeTrue«N+(Ql- 
Q2). 

The count observed by the circuit is only N counts. 

The "true" time difference is the number of internal 
master clocks in the interval times the "true*' period of the 
clock. 

DeltaTimeTrue equals M^trueperiod-M/freqGR. 

Thus have: 

tf^Ql-Ql) equals freqCS*M(freqGR (eq A.). 
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Re-arranging, 

Qf/M)%baseGR+driftGR) equals (Ql-Q 
^ 1) • {baseGR+dnfiGR)/M+baseCS+dHfiCS. 

{N/M)*baseGR-baseCS equals {VAf)*'{Q2-Q 

lybaseGR+driftCS- 

drmGR*(l/Af)*(-(!22-(?l)+A0. 

The term in the brackets can be replaced by manipulating 
eqA. 

15 

{N-¥(Ql-Q2))fM equals freqCSIfreqGR. 

Thus, 

Q^/Af)*baseGR-baseCS'- 
(l/M)*(Q2-Ql)*bascGR+ 
driflCS-driftGR*(rre?CS//«^/?)- 
25 drifiCS-driftGR *(fn:qCS(fhqGR)+{VAf)*iQ2-QiybaacGR. 

The scale factor in the first term on the left side of the 
equation (baseGR/M) is the approximation of the deltaTime 
based on the nominal GPS receiver frequency. Thus, this 

30 term is the estimate of the total cell-site frequency. By 
subtracting the nominal cell-site frequency from the total 
frequency estimate, the difference represents the difference 
in the frequency offset of the two clocks. The left hand side 
represents the manipulation of the measurement data that 

35 produces the frequency difference wanted. The error in the 
frequency estimate is due to the lost fractional cycles over 
the time interval. The error is reduced by either increasing 
the cell-site frequency (increasing baseCS which decreases 
Ql and Q2), or increasing the time interval (by increasing 

40 M). The scale factor on the GPS receiver drift term (fireqCS/ 
fireqGR) produces the translation of the GPS receiver clock 
frequency from its reference frequency into the reference 
firequency of the cell-site. 
The freqDiff measurement N is defined as shown below, 

45 

Af equals (NfAf)*baseGR-baseCS (units are counts at 
baseCSfnq) error equals {l/Af)'{Q2-Ql)*baseGR 

Applying this change of variables, 

50 

DriftCS equals AF+driftGR*(baseCS+driftCS)/{baseGR+drifiGR)- 
error. 

Expanding, 

55 DriftCS^ii-dHftGRIibaseGR+driftGR))" 

t^+dHftGR*baseCS/(baseGR+drtftGR)-errOK 
Expanding of the left, 

DriftCS*{(baseGR+dnftGR-driftGR)/{baseGR+driftGR))'' 
tsF-¥dHftGR*baseCS!{baseGR-^driftGR)-error. 

And reducing, 

DriftCS*baseGRJ{baseGR+driftGR)~ 
^F+driftGR*baseCS/{baseGR+dnftGR)-error. 
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Finally, solve for the cell-site drift and re -insert the error 
term, 

DriftCS equals AF*freqGR/baseGR-¥dnftGR*baseCS/baseGR+ 

Using known parameters, the final cell-site drift estimate, 

DriftCShat equals AJ^*freqGRhat/baseGR+tlnftGRhat*baseCS/ 
baseGR. 

Where AF equals (N/M)*baseGR-baseCS, 
In the cell-site drift estimate accuracy the difference 
between the "true" and measured ceU-site drift is 

errorDriftCSmerrorDrifiGR 'AF/baseGR+errorDriJtGR *baseCS/ 
bascGR+ 

{lfAf)^(Ql-QZ)*fKqGR. 

Inserting the definition of the AF, 

BrrorDriftCS^ 

cnoTDriRGR'(VbaseGIt)*(l^/M^baseGR-baseCS+baseCS)-i- 
{\iM)*(gUQiyfreqGR^ 
errorDriftGR*(l^haj:fiC?/i) • Q^/M*baseGR)+ 
{Vhf)*(Ql-Ql)yfnqGR« 
enorDriftOR*N/M+CQl-(22)/MVre<7(7/i= 
cnorf)rif tGR*N/M+((ei-G2)/A/) *frtqGR, 



Thus, 

ErrvrDrifiCS equals {{Ql-Q2)/M)*freqGR+i^IAf)'erroTDriftGR. 

The longer the integration time, the more the error in 
estimating the "true" time interval grows because of inte- 
grating the GPS drift error term errorDriftGR. The quanti- 
zation error is the sum of two uniform random variables. 

It is the period of the clock used to count the frequency 
difference that determines the worst case frequency error. 
The higher the frequency, the smaller the period and the 
quantization term. Alternatively, one can use a'longer obser- 
vation period. Such may be a problem if the clock operates 
in an environment of some observable temperature change. 
This may cause the clock to change frequency, and can lead 
to non- linear changes. The average frequency will be less 
accurate than the quantized estimate. 

If the external clock frequency is lower that the internal 
clock frequency, the driftCS estimate from the timeDifif 
measurement is used to estimate the cell-site frequency 
offset. However, if the external frequency is higher than the 
internal clock frequency, then the freqDiff observation is 
used instead. 

In cases where only frequency synchronization is 
possible, because a system IPPS does not exist, then the 
freqDiff counter will be used to estimate the extemal clock 
frequency as the timeDifif will not provide any information 
about the external clock when the extemal PPS is not 
available. 

For a device that communicates with a network infra- 
structure and a GPS receiver, both receivers will track the 
carrier frequency of their respective transmitters. For GPS, 
the receiver also tracks the phase of the pseudorandom code 
allowing a pseudorange measurement to be produced. The 
measurement is a time difference of arrival measurement, as 



the time of transmission is known, and the receiver records 
the time of arrival with respect to its local clock. 

The GPS satellite (SV) has a time and frequency error 
from "true" GPS time and the "true" desired reference 
frequency of the spacecraft clock. 

FreqSV equals baseSVfnq+drifiSV. 

TrnieSV equals timeGPS-^-biasSV 

10 The GPS receiver (GR) also has a time and frequency 
error from GPS time. 

FreqGR equals baseGRfreq+drifiGR. 
TimeGR equals UmeGPS+biasGR 

15 

A cell-site system (CSS) with a reference time for all its 
cell-sites uses cell-site^ystcm time or timcCSS. Such time 
reference is different from GPS time. It is also possible that 
the cell-sites are not synchronized to a common time. 

20 Embodiments of the present invention do not require that 
any of the cell-sites have the same time reference. 

"FreqCS" equals baseCSfreq+driftCS. "TimeCS" equals 
timeCSS+biasCSS (with respect to cell-site System time). It 
also equals timeCSG+biasCSG (with respect to GPS time), 

25 Each cell-site has an offeet from timeCSS that is the 
biasCSS. We can also reference time CSS to GPS time. The 
time error for the cell-site changes to biasCSG. The cell- 
phone also has a frequency oflfeetfrom its nominal frequency 
referred to as driftCP. 

30 

FreqCP-baseCPfrcq+drLftCP 
rimeCP-timeCPS+biasCPS-timeCPG+biasCPG. 

There are two cell-phone methods for tracking the cell- 

35 site frequency, the NCO method or the VCO method. In an 
NCO method embodiment of the present invention, the 
cell-phone reference frequency is held constant and the 
additional frequency required to track the cell-site is made in 
an additional oscillator or NCO. The circuit that produces 

40 this additional frequency will generally have a software 
counter that integrates all the frequency adjustments to 
produce a nimiber that represents the ofifeet from nominal 
frequency needed to track the cell-site frequency. Such 
frequency offset is the driven by three processes, (1) the 

45 cell-phone offeet from its nominal frequency, (2) the Dop- 
pler between the cell-phone and the oell-site caused by the 
projection of the cell-phone user dynamics along the line- 
of-sight between the cell-phone and the cell -site, and (3) the 
frequency ofket of the cell-site from its nominal frequency. 

50 The cell-phone outputs a message that contains a total 
frequency offset number and an associated time tag. Such 
number is used with the freqDiff and timeDiff observations 
to synchronize the cell-phone to the GPS receiver clock. 
In the VCO method embodiment of the present invention, 

55 a cell-site carrier is tracked by adjusting the reference 
frequency directly with a VCO. The voltage is changed to 
the circuit that produces the different mixdown frequencies. 
The NCO is not needed. If the adjusted reference clock is 
directly observed with the timeDiff and freqDiff 

60 observations, then there is no need for the cell-phone to send 
any additional information. Any cell-phone bias is replaced 
by a cell-site range and bias, and any cell-phone drift 
becomes the raugeRate and cell-site drift. 

The tempMeas observable is very similar to the freqDiff. 

65 The navigation digital-signal processor preferably includes 
an oscillator whose frequency is determined by an RC-time 
constant in which the resistor is actually a thermistor. The 
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tempMeas circuit counts the number of cycles of the oscil- 
lator in a predetermined interval and gets a measure of 
temperature from a thermistor calibration curve. 

Tlie GPS dock naturally drifts with temperature. Thus, 
during a calibration phase when the temperature of the 5 
complete navigation digital-signal processor based GPS 
receiver, observe that the count and the GPS receiver clock 
drift (driftGR2hat) is swept at the same time. A model is 
built that allows the GPS clock drift to be predicted from the 
tcmpMcas. 10 

Tlie aging of the clock and the crystal imperfections will 
cause deviations of the crystal frequency from the model. 
And the intrinsic variations in the crystal act to limit the 
accuracy of the estimator. 

The hardware to implement the tempMcas is simply a 15 
counter that starts on a GPS predetermined millisecond and 
ends after a predetermined number of milhseconds later. 

While sweeping over the industrial temperature range in 
the factory (and in the field whenever obtain a reliable GPS 
fix), observe tempMeas and the driftGR2hat and build the 20 
model (and update model in the field) of the driftGR2hat 
verses the counter value. A tempMeas can then be performed 
at any time and using this model, an estimate of the GPS 
receiver clock drift can be obtained. This estimate is called 
driftGR2HatSCX0. The acronym SCXO comes from 25 
Software-Compensated-Crystal-Oscillator. 

In the cell-phone, a tempMeas at startup to form 
driftGR2hat is made. Such estimate will be used during the 
time and frequency transfer (from the phone to the cell- 
phone) to reduce the required GPS search range and thus, 30 
enable a super-thin GPS client. 

The components of the dual mode system are a GPS 
satelUte, the GPS receiver, the cell-site, and the cell-phone. 

A GPS satellite is being tracked by a GPS receiver at the 
cell-site whose purpose is to observe the offset of the 35 
cell-site clock from GPS time and its frequency reference. 
This receiver and cell-site have known locations and are 
static. They are physically close enough to allow the physi- 
cal connection needed to make the timeDiff and freqDiff 
observations. 40 

The second GPS receiver is at the cell-phone location. 
The goal is to estimate the location (position and velocity) 
of the cell-phone. All possible measurement sources to do 
this. The normal GPS pseudorange and pseudorange-rate 
measurements are used to compute the position, time, and 45 
frequency of the GPS receiver at the cell site. Such allows 
the GPS receiver to synchronize to GPS time, and thus, by 
observing the relative ofiket between the cell-site and the 
GPS receiver time and frequency, can calculate the time and 
frequency offiset of the cell-site from GPS time. 50 

In cell-phone embodiments of the present invention, the 
inverse is done. The cell-phone is assumed to be tracking the 
cell-site frequency and time, Pseudorange and pseudorange- 
rate measurements can also be formed for the signal between 
the cell-phone and the cell-site. Because the cell-site posi- 55 
tion is known, the ceU-site can be used as a ranging source 
in the same way a GPS satellite is used. 

Since the cell-phone is tracking the cell-site, and the 
frequency required to track is known, the use of freqDiff 
allows two things. The pseudorange and pseudorange-rate 60 
between the cell-phone and the cell -site can be used as 
estimates of the cell-phone time and frequency error. Esti- 
mates can be made about the GPS receiver's time and 
frequency using the timeDiff and freqDiff observations, and 
greatly reduce the search range for the GPS satellites. It can 65 
also be used to improve acquisition time and/or acquisition 
sensitivity. The cell system aids the GPS to acquire satellites. 
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After at least one GPS satellite is acquired, the pseudorange 
and pseudorange-rate between the cell-phone and the cell- 
site are used along with the GPS pseudoranges and 
pseudorange-rates in the position calculation. Because of the 
timeDiff and freqDiff observations, the cell-site time and 
drift can be expressed in terms of the GPS bias and drift, and 
the number of unknowns is reduced to the normal unknown 
as in the stand-alone GPS case. 

FIG. 7 represents a dual-mode system 700 with GPS 
receiver tracking a satellite vehicle (SV) 702 and a cell- 
phone 704 tracking a cell-site 706. The clocks in the system 
with respect to the dual-mode system are, 

SV clock* (numSatelUtes), 

GPS receiver-1 clock in the cell-site, 

GPS receiver-2 clock in the cell-phone, 

Cell-site clock*(numCellSites), 

Cell-phone clock 

For example, if there are five GPS satellites, and five 
cell-sites, then the total number of clocks in the system, 
nmSatelLites+numCellSites+3=5+5+3=13. For each cell -site 
and GPS satellite, a pseudorange (PR) and pseudorange-rate 
measurement (PRR) can be formed. 

For a system with one GPS satellite 702 and one cell-site 
706, the calibration begins by making a passive connection 
between a cellnsite 706 and cell-phone 704. A cell-site clock 
signal and its 1-PPS can be observed by GPS rcceiver-1 706. 
These signals are cormcctcd to the freqDiff and timeDiff 
circuits of a navigation digital-signal processor respectively. 
The cell-phone clock signal and its 1-PPS can be observed 
by the GPS reo6iver-2 704. Such also has a navigation 
digital-signal processor chip and such signals are connected 
to the freqDiff and timeDiff circuits of navigation digital- 
signal processor respectively. 

FIG. 8 represents each clock with its own system time and 
a corresponding clock bias from this time. The sum of these 
two, e.g., the system time plus the bias, produce the refer- 
ence time used by the hardware for making observations in 
the receiver. The time-difference observations are taken 
between these time-references. The first one, TDl, is 
between GPS receiver-l and the cell-site. The second one 
TD2 is between GPS receiver-2 and the cell-phone. The 
clock pulse signals represent a timing signal leaving the 
transmitter at a known time and arriving at the receiver later 
in time according to the range. However, the time bias errors 
cause the observed phase measurement to be affected, e.g., 
by the bias error of both the transmitter and receiver. Such 
observations become the pseudoranges. 

TDl is the timeDiff measurement between the 1-PPS of 
the GPS receiver-l and the 1-PPS of the cell-site. Such 
produces the observation N12 and N21 in units of counts 
from the receiver-l and can be used to form timeDiff TDl 
in units of seconds. 

TD2 is the timeDiff measurement between the 1-PPS of 
the GPS receiver-l and the 1-PPS of the cell-site. Such 
produces the observation N12 and N21 in units of counts 
from the receiver-2 and can be used to form timeDiff TDl 
in units of seconds. 

From FIG. 8, it can be seen that, 

BiasORl+'roi«bia5CSS+A(CSS-GPS)-bia5CSO, (units are sec- 
onds) 

BiasGR2+TD2-biasCPS+A(CSS-GPS)-biaBCPG. 

The units of measure for the pseudorange measurements 
are preferably maintained in units of seconds rather meters, 
in order to simplify things in terms of the timeDiff and 
freqDiff. But the units in meters can be used too. 
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Two measurements are obtained from the timeDifif circuit 
each second. They arc Nil and Ni2, where I=timc marker, 
and "I" is the measurement from the external PPS and "2" 
is the measurement from the internal PPS. 

The phase measurement that was produced to make the 
local phase of the signal -generated in the cell-phone track 
the received phase of the signal transmitted by the cell-site 
can be accessed. 

From FIG. 8, the observed phase measurement is such 
that, 

BiasCS+rangeCP'BiasCF+^p,(umts are seconds) 

Re-arranging this yields, 

^p-rangeCP+BiasCS-BiasCP, 

TDl'^asCSG-BiasGRl, 

TD2^asCPG-BiasGR2 

The analogous measurement from the GPS receiver is 
such that, 

l){asSV+rang€GRl''btasGR2-H^R2, 
UasSV^rangeGRUUasGRl-^Rl, 

Re-arranging gives, 

i^RUrmseGRV^biasSV-hiasGRl, 
^Rl-rangeGRl^biasSV-biasGR2. 

The "true" positions can then be found. A three- 
dimensional earth-centered-earth-fixed (ECEF) Cartesian 
coordinate, e.g., WGS-84, is used. A cell-site position and 
velocity is considered to be known, 

/*asC5(m meters)-(rC5; yCS, zCS), 

\^lCS(m m/s)=(jVc/C5, yVklCS, zVelCS)«0, 

The GPS receiver-l at the cell-site has a known position 
that is near the cell-site position, but does not need to be 
identical to it. Thus, posGPSl«(xl, yl, zl), and velGPSl= 
(xVell, yVell, zVell)=0. 

The GPS satellite position is assumed to be known from 
the ephemeris model that is distributed around the system. 
The satellite velocity is calculated from the ephemeris. 
PosSVo(xSV, ySV, zSV), VelSVo(xVelSV, yVelSV, 
zVelSV). 

The GPS recciver-2 has a cell-phone and its loaction 
represents the main unknown in the system. All the other 
unknowns arc nuisance parameters in the system that need 
to be solved or eliminated. 

The cell-site and GPS receiver are defined to be at the 
same (unknown) location. So, posGPS2=posCP=(x, y, z), 
and vclGPS2=velCP-(xVel, yVel, zVcl). 

In the "true" range definitions, the coordinate frame is 
orthonormal. The positive distance between any two points 
is the magnitude of the vector between the two points. The 
range equations can be expresses in terms of seconds. 

Rang€Gri*k*((xSV-xiy+(ySV-ylf+{zSV-ziyy^. 

RangeGR2''k*{(xSV-xfi-(ySy-y)^+(zSV-f)'^. 

RangeCPmk*{{xCP-xf^{yCP-y)M^P'^f)^- 

Where, K=the inverse of the speed of light, e.g., 
2.99792458e8 meters/second. 
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FIG. 9 is a visualization of how the freqDiff measurement 
can be taken and used. A freqDififl measurement is obtained 
between GPS receivcr-1 reference and the cell-site refer- 
ence. The freqDifiE2 measurement is obtained between GPS 

5 receivcr-2 reference and the cell-phone reference. The same 
reference, e,g,, the nominal frequency used in the GPS 
receiver, is used for both the GPS receivers. Part of the 
strategy is to employ the same hardware solution that 
contains the navigation digital-signal processor chip at both 

10 the reference stations for DGPS and data, cell -sites, and 
client (ceU-phonc) solutions. Thus, the same frequency 
reference can be assumed, which is the nominal crystal 
frequency used to generate all the frequencies used to track 
the GPS satellite carrier. 

15 For the cell-system side, it is convenient to use a common 
frequency reference when calculating the frequency offset 
from systems with different nominal frequencies. However, 
the same modeling can be used in either case. All the 
nominal frequencies need to be scaled to one master fre- 

20 quency. Then the frequency offset is translated from its 
nominal frequency to the new master frequency by multi- 
plying the offset by the ratio, othcrRefFreq/nominalRefFrcq. 
An algebra can be performed with proper scaling. The result 
is converted back to the nominal reference with the inverse 

25 scaling. 

The range-rate measurements are not included in the 
freqDiff figure, the pseudoranges are in the timeDiff figure. 
The user motion and satelhte are usually expressed at the 
GPS carrier frequency (Ll-1575.42 MHz). The freqDiff 

30 observations are not affected by the user or satellite motion. 
Such applies to the NCO for the cell-phone frequency 
generation. In the VCO method embodiment of the present 
invention, the range and cell-site clock effects do appear in 
the cell-phone reference frequency. 

35 The raw measurement from the freqDiff circuit are the 
counts Nl from GPS receiver- 1 and N2 from GPS receiver- 
2. The freqDiff in units of cycles/sec at the cell-site reference 
frequency can be computed with these. BaseCS-baseCP can 
be assumed, 

40 

d^'-{NVM)*hascGR-baseCS (Hz at 
baseCSfreq) 

t^'lfaiM)*baseGR-baseCP (Hz at 

45 

baseCPfreq) 

where, 

NlofreqDiff in coimts from navigation digital-signal pro- 
cessorl 

N2«fr6qDiff in counts from navigation digital-signal pro- 
cessor2. 

An estimate of frequency from the cell-phone is required 
to track the cell-site carrier. The method embodiment of the 
55 present invention is used with an NCO. The cell-phone can 
output a number that can be scaled into units of cycles/sec 
at the base CP frequency. 

The measurement from the cell-phone carrier tracking 
loop, 

60 

^otCP^rangeRateCP+driftCS-drifiCP 

(Hz at baseCSfreq) 

In another embodiment, the VCO is steered to track the 
65 cell-site carrier and the effect on the cell-phone nominal 
frequency is observed directly in the timcDiff and freqDiff 
circuit. It is still possible that the numerical representation of 



03/04/2004, EAST Version: 1.4.1 



us 6,4 

25 

the frequency offset from nominal can be output. However, 
there is a different meaning of the frequency component due 
to the cell-phone dock that appears in frcqDiff. Now the 
freqDiff will observe the Doppler component of the cell- 
phone motion and the frequency offiset of the cell-site. 

If the transmitter has a positive frequency offset ^ then the 
receiver must increase its generated frequency to track this 
addition frequency offset from nominal. If the receiver has 
a positive frequency offset, then it must reduce its generated 
frequency to track the carrier Finally, if there is a positive 
Doppler between the transmitter and receiver, then the 
receiver must increase its frequency to track this effect. 

Tracking the carrier of the cell-site, the frequency pro- 
duced to allow tracking, 

FreqCPvco = 

freqCP(&x turn on) + {mngeRateCP ¥ drifiCS - driftsCF) - 
baseCPfreq-^ driftCP+ {rangeRateCP+ driftCS - drifisCP) = 
baseCPfreq + (rangeRareCP driftCS) = 

baseCPfreq 4- driftCPvco^ 

where, when tracking the cell-site, DriftCPvco« 
rangeRateCP+driftCS. 

The phase error discriminator that produces the updates to 
the VCO frequency must remove the drift of the cell-phone. 

Examining the effect on the freqDiff observation, the 
counter will observe the integer number of cycles of the 
frequency freqCPvco to track the cell-site, rather than just 
the cell-phone frequency freqCP. 

In other words, ClockCountTrue = 
freqCSvco*delta'limeTrueeNvco+(Ql-Q2). 

Thus, when equating the observation to the unknowns, 
define, 

FregDijJVCO = dr'^tCPvco- driftGR* (JreqCS } freqGR) = 

{rangeRateCP+ driftCS) - driftGR* [freqCS/ freqGR) 

The GPS receiver can make a measurement of the 
pseudorange-rate from the frequency offset it generates in an 
numerically controlled oscillator that provides an additional 
frequency offiset from the reference frequency (freqCS= 
baseGRfi-eq+driftGR). The measurement can be modeled as 
the derivative of the pseudorange, 

A "Dot" after ^ indicates the derivative. 

^otGRl'rangaRatcGRi+driftSV-driftGRl ; 
^0tGF2~remgeRateGRl+driftSV-driftGR2. 

Removing unknowns for the freqDiff and pseudorange - 
rates, driftSV is assumed to be known from the ephemeris 
data that is distributed around the system. The cell-site is 
known to be static, and rangeRateGRl can be computed 
using, 

R<mgeRateGm'(x\%lSV-xVt)*uX+{ 
yVc\S-y\ciy uY+(zVklSV-zVt)*uZ 

where, 

uX=(xSV-xl)/raQgeGRl 

uyo(ySV-yl)/rangeGRl 

uZ=(zSV-zl)AangeGRl and 
"true" user velocity vector=(xVel, yVel, zVel) 



3,030 Bl 

26 

The GPS receiver at the cell^ite can be assumed to have 
a "true" user velocity of zero. Thus, all the positions and 
velocities are assumed known. The GPS receiver 1L2 at the 
cell-site can compute its position, or the position can be a 
5 surveyed position. Secondly, the satelUte position and veloc- 
ity can be computed with the ephemeris made available to 
the cell-site GPS receiver. 

Because RangeRateGRl and driftSV are knovra, can use 
the observed pseudorange-rate can be used to estimate the 
10 GPS receiver drift. 

Thus, the estimate of the drift at GPS receiver-l can be 
computed. 

driftGRlhat^rangeRateGRl+driftSV-^DotGRl{unHs are 
Hz at baseGRfreq). 

If more satellites are being tracked, then the drift estimate 
can be formed from the velocity fix, or by averaging a set of 
the above equations if the position is known. Such will 
20 improve the accuracy of the estimate if the position accuracy 
is high. 

Using the freqDiff and the estimated GPS receiver-l clock 
drift, the estimate of the cell-site clock drift can be formed. 

25 DrifiCShat-AFl *freqGR\hca}hascGR^ftGRl 

hat'baseCS/baseGR (units are HZ at baseCSfivq) 

where, 

AFl=(Nl/M)*baseGR-baseCS 

freqGRhatl=baseGRhat+driftGRlhat (from the GPS 
receiver-l) 

Nl "freqDiff between cell-site and GPS receiver-l. 
Remove the cell-phone clock drift to compute velocity 
35 with method embodiment of the present invention. 

When the cell-phone is tracking frequency with an NCO, 
freqDiff can be used to express the cell-site drift in terms of 
the GPS receiver drift. This removes the cell-phone clock 
from the measured cell-phone carrier frequency to track the 
40 cell-site carrier. 
Remember: 

= (yV2/ M) * baseGRfreq - baseCPfreq = 
45 driftCP - driftCRl * {baseCSfreq / baseGRfreq). 

Thus, 

DnftCPhat^hFl-^driftGRl^ihaseCSfnqlbaseGRfnq) 

50 

The carrier measurement can be used as an estimate of 
driftCP: 

^0tCF'^ngeRateCP+driflsCShat-(AF2+driftGRl'( 
ba8eCSfreq/baseGRfieq))(Hz at baseCSfreq). 

Thus, this is an expression for the pseudorange-rate where 
the unknown are in terms of velocity and clock bias of the 
GPS receiver. 

60 Remove the cell-phone clock drift to predict the drift of 
the GPS receiver in the cell-phone. 

Before fixing, this procedure can be used to generate an 
estimate of the clock drift of the GPS receiver clock in the 
cell-phone. 

65 An estimate of the drift can be formed by assuming that 
the cell-phone is static, so that the rangeRateCP term can be 
assumed to be zero. The accuracy of the estimate is driven 
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by the accuracy of this assumption. However, for most low The "truth" model for the observatioo, 

cost clocks, the frequency uncertainty caused by tempera- 
ture and aging od the crystal are much larger than the user FreqDiffVC0'^((NvcxH-(Ql-Q2)IAf)*freqGR (in 
dynamics between the cell-phone and cell-site, and the 
estimate provides a huge information gain that can be used 5 

to reduce search time and or increase sensitivity^ ^^^^^ „^ ^^^^^ measurement. 

For example, suppose the user is moving at 100 mph. We 
first convert the user velocity to meters/second: FreqDiffVCO Cobservcd)-Arvco (in integer counts at CP 

(100 miies;/hour)'(l -hour/3 600 seconds) •(5280 feet/1 -mile) jq &eq). 

•(12 inches/1 foot)*(D.0254 metei8/i-iach)-44.704 metcra/second Wc then form the freqDiff estimate in units of Hertz. 

Now, for this example, assume the nominal frequency of ^vco~{Nvco/M)'baseGiifnq~baseCSfKq 
the GPS receiver at the cell-phone clock if 27.456 MHz. 

15 (counte at CP freq). 

The wavelength of the. 



XCP=C/f=speed-of-light/27.456e6cycle/scc-10.9 
meters/ cycle. 



Our model is then, 

lJ^vco<'drtfiCSvco-driftGR2*{baseCSfreq/baseGRfnq) 



(in Hertz at CP (Teq)''(rangeRateCP+djiftCSy 

The frequency Error of drift CP assuming the user is static dnfiGRl*{baseCSfreqibaseGRfreq). 
when in fact the user is moving at 100 mph, , 

Inserting the estimate of the ceU-sitc drift, a measurement 
freqError (Hz at nominal GR freq)^(AAJQ4 mctei8/second)'(l- equation that can bc used to solve for the user velocity and 

cycle of 25 GPS receiver clock drift is made. 

baseCPfreqnO.9 meterB)-4.1 Hz. ^Ylvco«rangeRateCP^drifiCShat-driflGR2*{baseCSfreqlhaseGR- 

fnq) 



The freqError in Parts-Pcr-Million (PPM), 

fre^mwCPPM)=4.1/27.456=0.15 PFM. 



(in Heitz at CP freq). 

30 



Remove the cell-phone clock drift to predict the drift of 

The effect on the velocity solution is that this frequency ops receiver in the cell-phone 

error is multiplied up to LI whose waveleDglh=0.19 meters/ ^he user velocity is zero. Such allows an estimate of the 

^y^* clock drift of the GPS receiver in the cell-phone. 

Assume that the rangeRateCP is zero, 

Frequency Error at Ll-(44. 704 mctcis/sccond)*(l -cycle of ^ -r^^^, v ^ 

Ll/0.19 meters)»235 Hz. DnfiGRlhat^{dHftCShat-AFlvcoYbastGRfTtqlbaseCSfrcq 

Hie PPM remains the same at any reference frequency. 

The freqError in Parts-Per-MilUon (PPM), 235 Hz/1575 .42o ^ ™ • , , ^ . ^ .x. c 

0 15 PPM estmaate can be used to reduce the frequency search 

■ The frequency error of the crystal may be on the oirier of '"'^'^'^^ ^"'5"'"'=^ 

0.5 to 1 PPM, and the assumption of zero-user-velocity still '^"^ ''^ performed with this observation. 

provides a big information gain in estimating the GPS ^° ^^"^"^^ measurement for 

receiver frequency drift. The information gain for a 1 -PPM ^ ^ ^- . c.^ , ■ ^ 

• - ,r^\. i- ^ n m • . ^ . method embodiment of the present invention (using the NCO 
uncertamty is 1/0.15=6.7. This means the frequency search 

can be reduced by a factor of 6.69. Such provides a major b the cell-phone) is <t>DotCP-rangeRateCP+driftaCShat 
reduction in search time. 

Assuming the user dynamics are zero, first make an -(AF2+driftGR2*(bascCS&cq/bascORfreq)). 
estimate of the cell-phone frequency drift, 

Also the measurement for another embodiment of the 

DriftCFhat^rangeRateCF+driftcshat-^otcp present invention (using the VCO in the cell-phone) is 

(Hz at baseCSfreq) ^U^vco^rangeRateCP-HiHftCShat-driftGR 

And then, using the freqDiff and assuming the l-{bastCSfrtqiba5eGRjrtq\ 
rangeRateCP=0, estimate the GPS clock drift as, 



DnftGR2hat'{driftCPhat-AFZ)*(baseGRfreq/baseCSfreq) 

(Hz at bas€GRfreq)~{dnftCShat-i^DotCP-AF2)*{baseGRfreq/ gQ 
baseCSfreq) 



These equations are generalized for simplicity as, 

PrrCPmC3*rangeRateCP-driftGR2*Ct+C2 



where for a method embodiment of the present invention, 

(units are Hz at baseGRfrtq), PrrCP-(t)DotCP 

Remove the cell-phone clock drift to compute velocity. In Cl=^aseCSfreq/b^eGRfreq) 

the case that the cell-phone is using the VCO to track the 65 C2=dnftsCShat-AF2 

cell-site frequency and the freqDiff can observe the effect on C3= 1/(XCP) 

the baseCPfreq. }.CP=C/baseCSfreq. 
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For another method embodiment of the present invention, 
PrrCP''£iF2vco 

Cl=(bas6CSfrcqybas6GRfrec|) 

C2=driftsCShat 

C3=l/(XCell-site) 

X-C/baseCPfreq. 
Combine the ccU-phone and GPS pseudorangc rate mea- 
surements to estimate the user velocity and, the GPS 
receiver drift. 

The cell-phone can track a set of J cell-cites and form the 
j-th pseudorange-rate observations as: 

PrT€P(j)«C3*r(mgeRateCF(jydnfiGR2*Cii-C2{fiz at baseCP- 
freq) 

K GP*S satellites form the k-th pseudorange-rate obser- 
vations as: 

PrG]a(k)~CS*rangeRateGR2(k)'driftGR2*CMS(i{z at baseGR- 

where: 
C4-1 

CS=driftSVhat (from the satellite ephemeris) 

C6-1/(>.GR) 

()^-GR)=C/baseGRfreq 

In order to solve for velocity, first perform the position 
solution, and then, predict the nominal rangeRateCP(j) as 
follows, 

RangeRateCPhat(f)-xVelCS{k)*u?f(k)+yVelS 
{kyuYik)+zVelCSik)*uZ(k), 

where, 

uXCS(j)-(xCSO)-x)/rangeGR 2 
uYCS(j)=(yCS(j)-y)/rangeGR 2 
uZCSO)-(zCSG)-z)/rangeGR 2. 

However, for the cell-site, assume it is static so, xVelCS 
(k)-yVelCS(k)-zVelCS(k)-0. 

Similarly, can predict the nominal rangeRateGR2 (k) as 
follows, 

RangeRateGKlhat{ky 
(xVelSV{k)-x\%t) *uX(k)-i-(yVeiS(k)~ 
y Vcl) *uY{k)-i-izVelSV{k)-zVkr)*uZ{k), 

where, 

uXSV(k)-(xSV{k)-x)/rangeGR2 
uYSV(k)=(ySV(k)-y)/rangeGR2 
uZSV(k)=(zSV(k)-z)/rangeGR2. 
Making these substitutions, 

UnPrrCP{j) = PrrCP{j) -C3* range RateCPhai{j) = 

drifiGRlMCJ A-C2 
UnPrrGR2{k) = PrrGR2{k) - C3 *rangeRareCR2haKk) = 
C6 * (-uXSV{k)*xVel - uYSV{k) * yVel - wZSV(fc) * zVeC) - 
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-continued 

driftGR2*C4+C5, 

5 These are the equations in terms of only the user velocity 
and the GPS receiver clock drift. Thus, as long as four 
measurement from either multiple GPS satellites or multiple 
cell-sites are found, the velocity and drift using conventional 
methods (like least-square or a Kalman filter) can be solved. 

To removing unknowns for the timeDifif and 
pscudorangcs, the timcDiflf measurements must be used. 

The rangeGRl is considered to be a known quantity since 
its components are known. 

The satellite position and clock bias is considered known 
since it can be computed from the GPS ephemeris that can 
be circulated around the system by the server. 

Assume that GPS rcceiver-l is tracking enough satellites 
to compute at least its time solution. If its position is known, 
then only one satellite is needed to do this. 

In timeDiff and pseudorange measurements, the pseudo- 
20 range are formed by the cell-phone tracking the cell-site, 
(j)CP=rangeCP+BiasCS-BiasCP. 

The pseudorange formed inside GPS rcceiver-l at the 
cell-site, (|)GRl-rangeGRl+biasSV-biasGRl. 

The pseudorange formed inside GPS receiver-2 at the 
25 cell-site, (t)GR2orangeoGR2+biasSV-biasGR2. 

With timeDiff observations between the GPS receiver-1 
and the cell-site (Nil and Ni2), biasCSGoBiasGRl+TDl. 

With timeDiff observations between the GPS rcccivcr-2 
and the cell-phone (Nil and Ni2), biasCPG-BiasGR2+TD2. 

To solve for the clock bias of the GPS receiver as the 
cell-site, the rangeGRl is known, and biasS V is known, and 
(JjGRI is measured in GPS receiver-1, can estimate the clock 
bias for GPS receiver-1 as foUows, 

btajGKl/wfcestimate of (biasGRl)'rangeCRl+ 
bLasSV-40Rl. 

To solve for the clock bias of the cell-site from GPS time, 
use the timeDiff measurement (TDl) between GPS 
receiver-1 and the cell-site. Estimate the bias between the 
^ cell-site and GPS time. 

Either of the timeDiff measurements (Nil or N12) can be 
used. 

TDUmVfTr^gGRihat 
TDl^NU/freqGRVuU 

where 

FreqGRhat=baseGR£req+driftGRlhat 
Assume that biasGRlhat and driftGRlhat are known and 
arc obtained directly from GPS receiver-1 at the cell-site. 
Now, using the definition of the timeDiff, solve for the 
50 cell-site drift as follows, 
using Nil, then 

BiasCSihat^TDl+blasGRlhat {units art seconds), 

55 using N12, then 

BiasCS\hat''l-TDl^biasGR\hat 
(tmits are seconda) 

Remove the ceU-phone clock bias to compute position. 
60 Assume the cell-phone produces the pseudorange mea- 
surement between the cell-phone and the cell-site. 

The pseudorange formed by the cell-phone tracking the 
cell-site is ((iCP^rangeCP+BiasCS-BiasCR 
Apply the estimate of the cell-site clock bias and use the 
65 timeDiff between the cell-phone and GPS receiver-2. Such 
change of variables allows elimination of the cell-phone 
clock and cell^te clock, and to express the cell-phone 
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derived pseudoranges in terms of the ccU-phone location and 
GPS receivcr-2's clock. 

The timeDiff in the cell-phone has the form, 

7Di~NiVfK(PRlhat 

TDl-Nll/freqGJWtat 

where, 

FreqGRhat«baseGRfreq+driftGR2hat 

Forming the timeDiff without an estimate of the driftGR 
in the cell-phone, do not know have an estimate of the GPS 
receiver-2 drift until a position fix. Thus, re-formulate the 
model to obtain a model that represents the explicite effect 
of this paramter. 

Normalize the observed timeDiff count with baseGRfreq 
rather than FreqGRhat. 

In this case, 

TDJ = NJJ/freqGR = 

N}I*{ij {baseGRfreq + driftGR) - 1 j baseGRfreq + 
1 1 baseGRfreq) = 
NJ J * {{baseGRfreq - baseGRfreq - dr0GR) { 
{freqCR* baseGRfreq) -t- 1 (baseGR) = 
Nll*{i{ baseGRfreq - drifiGRj {JreqGR * baseGl^req)) = 

Nl J /baseGRfreq - NJ I mdriftCR/ {freqGR* baseGRfreq) 

The measurement is not a linear fimction of driftGR 
because it represents up in both the numerator and denomi- 
nator (inside freqGR) of the second term. While this can be 
handled, it is inconvenient. Thus, perform a simple approxi- 
mation that linearizes the model. Replace the freqGR in the 
denominator with simply baseGRfreq. 

Linearize the unknown GPS receiver drift to make the 
measurement linear in the variable driftGR. Remove the 
variable driftGR from the denominator. 

Thus, examine the approximation of 

dHftGRj{baseGEfreq*freqGK)= 
di\&.QV^J(baseGKfrcq*baseGEfrtq), 

For again that driftGRo20PPM. 

Solving for the "true" value, = 

349/(27.456^6 * 27.456549tftf) = 7.28264497& - 13. 

If approximate the denominator as baseGR* baseGR, then 
the ratio, 

then the ratio, = 549/(27.456*6 •27.456e6) = 7.282790418* - 13. 

The difference=1.456220e-17. 

Multiply by the max Nl=baseGR, the error is 3.998197e- 
10=0.11 meters. 

Thus, there is still a linear error model. If this technique 
didn't work, then the model would not be linear since the 
driftGR term represents up in the denominator. 

To use the cell-phone pseudorange for positioning before 
have an estimate of driftGR2 is found, then include the drift 
term in the pseudorange. Normally, the pseudorange is only 
related to the position and clodc bias. Thus, this formulation 
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has an extra unknown. However, the pseudorange rates are 
available, all the unknowns can be solved together. 

To use the tempMeas to estimate driftGR2, suppose that 
driftGRo20 PPM. Suppose that Nil was almost exactly 
5 equal to bascGR. That means the cell-site is just to the left 
of the GPSl-PPS. The additional term proportional to 
driftGR is driftGR/frcqGR. 

If the drift were 20 PPM, then driftGRa20*baseGR/le6o 
549 Hz and driftGR2/baseGR+driftGR=0.2 seconds. 
10 The tempMeas observable are used. Using the SCXO 
calibration, the driftGR to an accuracy of+/0.5 PPM can be 
estimated. 

The driflGR=0.5*baseGR/le6 13.7 Hz. Such leads to a 
time enror of 13.7/(27,456,01 4)=5 .067=50 microseconds. 

15 Such is about 50 chips of the GPS pseudorandom code that 
repeats every 1023 chips. 

The tempMeas will provide a good enough approximatioD 
for using the ceU-phone pseudorange to predict the GPS 
receiver clock bias. Such can be used for prcpositioning. 

20 To work timeDiff measurement, we redefine the timeDiff 
for the cell-phone side as follows, Recall the measurement 
model, 

TDl^mNlVbaseGRfr^q-mX *driftGn/(freqGR* baseGR) 

^ Define, 

TDU^Nll/baseGI^q 

Thus, 

30 

TDU'TDl+}m*driftGRi(fKqGR*baseGR) 

Now, inserting the measurement model for the timeDiff at 
the cell-phone (don't include the quantization terms, so the 
35 parameters are the estimates and not the "true" parameters), 

TDX'biasCPhat-biasGKlhca. 

Thus, 

^ TDla^blasCPhat-blasGRlhat+Nll *driftGE2l 

{baseGRfreq *baseGRfreq) 

Using this form, the additional term proportional to 
driftGR2 will account for the error in using baseGRfreq 
45 rather than freqGR to normalize. 

Re -arrange the measurement model to solve for bias- 
CPhat. 

BiasCPhat=NlVbascGR^biasGRlhat-Nn * 
50 dz\5iG^{baseGRfreq*bascGRfniq), 

The above formula is inserted into the cell-phone pseu- 
dorange so that all the clock terms are in terms of GPS 
rec6iver-2's clock. 

55 

^CPmrangeCP-i-BiasCShat- 

{Nn/baseGR+biasGRlhat-mi*dHftGR/{baseGRfreq*baseG}^req) 

Define, 

60 

PtCPO^-^CP' 

m*rangeCP{f^biasGRl*DUDl (units 
of seconds) 

65 

Where, 
Dlol 
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10 



mVbaseGR~Nll*(driftCShat-<^otCP-An)/ 
{baseCSfrvq*baseGRfreq). 

The accuracy of the estimate is domioated by the assump- 
tion that the distance to the cell site is 1000 m. Even at a 
distance of 2000 m, the accuracy would be about 1000/c- 
3.3e-l-3.3 microsecond. The time-uncertainty has been 
reduced from 1 msec to 3.3 microsecond. That is an infor- 
mation gain of almost a factor of 300. Such is the kind of 
information gain that will allow a super-thin GPS client. 

Remove the ceU-phone clock bias to compute position. 

Derive the measurement models when the cell-phone uses 
the VCO method embodiment of the present invention to 
track the carrier. The frequency is adjusted to match the 
phase of the incoming signal. Tlie meaning of the timeDifif 
changes since the PPS associated with the cell-phone refer- 
ence also moves as the frequency is adjusted. 

The timeDiff no longer observes the cell-phone bias, but 
rather, it observes the cell-site bias and the range between 
the cell-site and the cell-phone. For example, if the cell-site 
clock bias is positive, then the cell-site must also increase its 
phase to be in alignment. To compensate for the travel time 
of the signal from the cell-site to the cell-phone, the cell- 
phone must further delay its signal in order to track the 
incoming signal. A phase advance is a positive increase in 
the phase, 

Tlie cell-phone clock bias is replaced by the rangeCP and 
the cell-site bias. 

The "true'* timeDiff models is: 



15 



D2obiasCShat-Nll/baseGR + Nll*driftGR/ 
(bascGRfreq*baseGRfrcq) (seconds) 

D3-1/C (units of seconds/meter) 

Remove the cell-phone clock bias to predict the bias of the 
GPS receiver in the cell-phone with a method embodiment 
of the present invention. 

To estimate the GPS receiver-2 bias directly, simply 
re-arrange the above formula to solve for it. 

biasGRlhat^rangeCP+BiasCShat-ipcpi- 

(NlVbaseGR-Nii*drifiGR/(bas€GRfirq*baseGRfreq)). 

This equation gives a method to estimate the GPS receiver 
clock bias prior to tracking any GPS signals. For GPS, the 
maximum time-search window is one -millisecond. Thus, in 
order to reduce the search range, a clock bias estimate with 
an error that is significantly less than one millisecond is 
needed. 

Suppose that the cell-phone is only one kilo-meter (1000 
m) from the cell-site. Such is possible in Japan where the 
ceU-site density is quite high. 

Also suppose using the freqDiff to estimate driftGR2. 

Recall, 



AF2^driftCP-dHftGR2 * {bastCSfreqIbaseGRfreq) 25 
DrifiGRZhat«{driftCShat-^otCP-AF2y{baseGRfKg/baseCSfreq) 



In this equation assume that the user velocity is zero. . 

Use the tempmeas to estimate driftGRscxo, that is, the 
estimate of the GPS receiver-2' s clock drift using the 
tempMeas observation and the calibration model to estimate 
clock frequency from the tempMeas. 

Assume the range from the cell-site is 1000, and use the 
driftGR estimate, 

biasGR2hatm(lO0/c)+biasCShat-^p^ 



where, 

BiasCPvcoorangeCP+biasCS 

We have the same issue for the actual measurement, i.e, 
to normalize the timeDiff by the correct frequency. Normal- 
ize baseGRfreq and then account for the driftGR2 tenm 
explicitely. 

Define, 



TDla^'NlVbaseGRSreq. 



Thus, 



TDU^TDl +Nn*driftGR 
{{baseGRfreq*baseGRfreq). 



Now, inserting our model for the "true" TDl, and insert- 
ing the estimated parameters, (not including the quantization 
terms in our measurement model), 



TD U-rangeCP-^biasCShat-biasGRl +Ntl*dnftGR 
l{baseGRfreq*baseGRfreq). 



30 



35 



40 



SO 



55 



Define, 

PrCP{))mNnjbascGR^D^*rangeCP{fy- 

biasGR2*DUD2 

(units of seconds) 

Where, 
Dl-1 

D2obaseCShat-Nll/baseGR + Nll* driftGR/ 
(bascGRfr6q*bas6GRfreq) (seconds) 

D3-1/C (units of seconds/meter) 

Have an expression for the cell-phone pseudorange where 
the unknowns are the user position, and the GPS receiver 
clock errors. Even though the GPS receiver drift appears, the 
equations can still be solved. 

Remove the cell-phone clock bias to compute position. 

For time transfer for aiding the GPS acquisition, simply 
re-arrange the above equation to solve for biasGR. Once 
again, have two sources to estimate the GPS receiver clock 
drift. Either use the tempMeas observation or use the 
freqDiff and assume the user velocity is zero. 

Biaj:GR2=rangeCP+biasCShat-mVbaseGRfnq+ 
mVdriftGR2f(baserGRfnq*baseGRfTTq). 

Use the data from the tempMeas observation, 

UasGKlha t-rangeCP+biasCShat-N 

lUbaseG]^Kq+Nll*driftGR2SCXO/ibaserGRfreq*baseGRfn:q). 

The combination of the freqDiff at the cell-phone assumes 
the cell-phone rangeRate is zero, 

BriftGRlha t'{dnftCSha t-AF2 vco) * baseGRfreq/baseCSfmq 

(in Hertz at GR freq). 

Inserting this, 

biasGRI^rangeCPA-biasCShat-m l/baseGRfreq*^N 
U*{driftCShat~^vco)/(baserCSfreq*baseGRfniq) 



TD I mBiasCPvco- biasGRl , 



(units are seconds) 
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To solving for PositionA^elocity in the Dual mode system. 
Combine the cell-phone pscudoranges with the GPS pseu- 
do ranges. 

Define for a method embodiment of the present invention 
(NCO method embodiment of the present invention), 

PrCP(f)-fpCJ'''D3 *rangeCP<J}- 

biasGR2*DUDl (units of seconds). 

Where, 
Dlol 

D2«biasCShat-Nll/baseGR + Nll*driftGR/ 

(baseGRfreq*baseGRfreq) (seconds) 
D3«l/C (units of seconds/meter). 
Or define for method embodiment of the present invention 
(VCO method embodiment of the present invention), 

PT€P{j)-m.HbaseGR'DZ*rangeCP{j)- 

biasGR2*Dl+D2 (units of seconds). 

Where, 
DUX 

D2 = baseCShat-Nll/baseGR + Nll *driftGR/ 

(baseGRfreq*bas6GRfreq) (seconds) 
D3=l/C (units of seconds/meter). 
And define the GPS observations as, 

PrORZ {kyD6*rangcSV{k)-DA'dnftGR 

2+Z>S. 

Where, 
D4=l 

D5=biasSV 
D60I/C. 

As with the pseudorange -rates, linearize them with 
respect to position by subtracting the nominal range. 

Define the delta vector as the error between the "true" 
position vector and the starting estimate of the position. In 
the estimation process, compute the delta vector and then 
add it to the starting position to get the position estimate. 
Thus, assume that when forming the predicted quantities 
where the user position is needed, use xNominal. 

Define the error in oxir nominal estimate (xNominal) as 
delX. 

delX-x-xNominal 

delY= y -yNom inal 

delZ=z-zNominal 

Assume having measurement to J cell -sites and K satel- 
lites. The linearized measurement for the j-th cell-site and 
the k-th GPS satellite are, 

UaPrCnj) = PrCP{}) - D3*rangeCPhaKj) = 

D3 * i-uXCSij) * delX - uYCSij) * delY - uZCS{j) * delZ) - 

bi(uGR2*Dl ^D2 

UnPrGR2{k) = PrGR{k) - C3 * rangeGRhaj{J) = 

D6 m i-uXSVik) *delX - uYSVik) •delY- u7SV{k) * delZ^ - 

blasGR2*D4^D5 

where the uXCS and uXSV arc the same as defined in the 
pseudorange- rate section. 

The Dual Mode system has derived the solution equations 
that relate all observable in the system to the user position/ 
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velocity and clock bias and drift of the GPS receiver in the 
cell-phone. All other system clocks have been cahbrated out 
of the system. 

Standard solution techniques can then be used to solve for 

5 the eight unknowns. 

To get a solution, many observations have unknowns. 
That means that four observation that have the user position 
are needed. The derivative measurements arc also needed. 
Estimate the GPS receiver velocity and clock drift. 

10 The cell-phone measurements along with the timcDiff, 
fi-cqDiff, and tcmpMeas obscrvablcs can be used to improve 
the acquisition of the GPS satellites by reducing the search 
range for those signals. Estimate the GPS receiver clock bias 
and drift with respect to the cell-phone and the calibration of 

15 the cell-sites. 

The system is a truly distributed processing system. 
Depending on the client type, all the calibrations and cal- 
culation could cither be done in the cell-phone, in the 
network server, or spread between the two systems. 

20 However, the server is the conduit that sends the calibration 
numbers for the cell-sites so that the data can be used for 
both time and frequency transfer that improves GPS acqui- 
sition times, and for Dual mode operation where use the 
range and range-rates from both systems to compute the 

25 cell-phone position. 

A communications means between the cell-site, the 
server, the cell-phone, and the GPS receivers at both the 
cell-site and the cell-phone is assumed. However, the details 
of the specific methods used are not critical to the equations 

30 presented here. 

GPS receiver-1 is at the cell site. It tracks GPS continu- 
ously and estimates its position, clock bias and clock drift. 
It performs a timeDiff and freqDiff between the GPS 
receiver and the Cell-site clock. It sends the measurements 

35 and its solutions back to the network server. Server com- 
putes the cell-site clock bias and drift. It also send other GPS 
data to the network server such as the 50 bps Navigation data 
message. 

The network server accepts the data from the cell-site and 

40 builds a table of the clock bias and drift of each cell-site. It 
also can generate a cell-site position from the GPS receiver 
position. The network server builds a Navigation Data 
database that is distributed throughout the system as a source 
of data for computing the satellite position, velocity, clock 

45 bias and drift. 

GPS receiver-2 is at the cell-phone to compute the posi- 
tion of this receiver. At turn on, don't know the GPS receiver 
clock bias, drift, or position. It performs a timeDiff, freqDiff, 
and tempMeas. The cell-phone receiver tracks a cell-site (at 

50 least one). By tracking the phase and frequency of the 
carrier, it can also form pscudoranges and pseudorange-rates 
between the cell-site. Using these measurements, plus the 
cell-site data that includes the cell-site clock bias and drift, 
and the cell-site position, and using the timeDiff, freqDiff, 

55 and tempMeas observations, can perform time and fre- 
quency transfer between the cell-phone and the GPS 
receiver. After acquitting one GPS satellite, can also use the 
observations to compute the dual-mode positionA^elocity, 
The time and frequency transfer allows us to improve the 

60 acquisition of the GPS satellites by reducing the search 
range. Such improves acquisition time and sensitivity. Such 
allows the super-thin-client. 

After at least one GPS satellite is acquired, the range and 
range-rate measurements from both systems can be com- 

65 bined to compute the cell-phone position. Such should allow 
a positioning capability when fewer than normal of the 
measurements from each system arc available. For example, 
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a three-dimensional fix can be made with only 1 GPS 
satellites, and 3 cell-sites, or 2GPS satellites and 2 cell-sites, 
or 3 GPS satellites and only one cell-site. 

Actually, one can compute the cell-phone position with- 
out any GPS satellite tracked at the cell-phone. Because the 
ceU-sites do significant relative vertical displacement, the 
three-dimensional solution is very noisy. However, a two- 
dimensional solution is possible if altitude of the cell-phone 
is available. (Remember, the GPS satellites are still needed 
at the cell-sites for synchronizing the cell-sites.) A three- 
dimensional-CDMA-only is also possible if the cell-sites 
have good vertical separation. 

GPS-navigation almanac and ephemeris are acquired and 
reduced by the webserver 106 with its own reference-station 
GPS -receiver to supply simple polynomials that respectively 
represent each sateUite's position and velocity. Such poly- 
nomials will, of course, go stale over time, but periodic 
contact between the ceU-site 102 and webserver 106 via the 
Internet are critically depended upon to provide fresh- 
enough orbit data. The GPS receiver 112 is therefore 
relieved of the storage-intensive requirement to hold current 
ephemeris and almanac data in local RAM/ROM, and 
relieved of the processor-intensive job of locally computing 
double-precision floating-point math the satellite-orbit posi- 
tions for every position fix. 

In a business-model embodiment of the present invention, 
the webserver 106 is owned and operated by an independent 
service provider. Per-use or subscription fees are charged to 
users who operate the cell-phone 104 and/or cell-site 102. 
The cell-phone 104 is alternatively further bundled in an 
intellectual property (IP) package for semiconductor inte- 
grated circuit device designs of application-specific network 
clients. Such IP is sold for a license fee and hardware-device 
cost to original equipment manufacturers (OEM's) of net- 
work clients. 

In operation, the GPS receiver 112 must simultaneously 
search for GPS-microwave signals in two domains, e.g., 
frequency and code-phase. The local oscillator and Doppler 
shift caused by satellite vehicle relative velocity create 
carrier frequency uncertainties that are resolved. The instan- 
taneous GPS-satellite pseudo-random number (pseudo- 
random number) code phase is another unknown. Received 
signals that are above the "noise floor*' are relatively easy 
and quick to search. But weak signals, as exist inside 
buildings, are buried in as much as twenty decibels of noise. 
Each visit to a frequency/code -phase bin must dwell there 
long enough to "beat down" the noise floor with processing 
gains provided by code correlators. Weak signals also 
require that the search bins have finer steps between bin 
frequencies and bin code-phase, e.g., due to aliasing. So 
more bins are needed to be searched and each bin needs 
more processing dwell, all of which increases search time 
exponentially. 

Hie webserver 106 publishes real-time navigation data to 
the Internet. Other services currently exist that put up partial 
or non-real-time navigation and initialization data on the 
Internet, and these can be used partially in place of, or to 
check the validity of real-time navigation data being sent to 
the GPS receiver 112. A primary and critical role job for 
webserver 106 is the ofif-loading of navigation- 
computational chores from the GPS receiver 112. Tlie 
degree to which such off-loading can be carried depends on 
the regularity of communication contact that can be realized 
between the GPS receiver 112 and the webserver 106. The 
actual measurements obtained by the GPS receiver 112 are 
forwarded to the webserver 106. 

Alternative embodiments of the present invention incor- 
porate webservers 106 that do database processing for 
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higher-level abstractions and purposes. For example, the 
server platform includes a health and quality monitor for 
checking said static observations and preventing an inclu- 
sion of incorrect information in said database of measure - 

5 ment errors and satellite data messages. On another front, 
webserver 106 could collect position solution information 
over time for the pattern of locations that GPS receiver 112 
visits over a period of hours, days, weeks, months, years, etc. 
Such information can then be processed to estimate where 

10 the user of GPS receiver 112 is, where the user has been, or 
where the user is likely to be in the future. Such information 
is useful to dispatch, time-clock, house- arrest, deployment, 
inventory, asset management, military, security, and other 
kinds of applications. Position information can be inter- 

15 preted to intelligently guess at what the user is doing at any 
one moment. For example, if the user's position matched 
that of a local grocery store, the user could be assumed to be 
shopping. Logs can also be generated by the user at the 
webserver 106, e.g., by city maintenance department oollo- 

20 eating with and "marking" in the electronic database the 
permanent location of a fire hydrants, power transformers, 
roadways, etc. 

The presence on the Internet of information about a user's 
position can be used by advertisers and marketers to direct 

25 contextual messages to the user in real-time. Such data is 
sold in real-time in a business model embodiment of the 
present invention. 

In general, embodiments of the present invention use 
Internet webservers to process data collected by remote GPS 

30 receivers in such ways that signal sensitivity is increased and 
search times are shortened compared to conventional receiv- 
ers and methods. 

Recent test data indicates the a priori worst-case modeling 
error for a typical crystal frequency is around ±0.5 PPM 

35 (±787 Hz) with temperature modeling. To be conservative, 
and leave some uncertainty for user position error effect on 
the frequency, a frequency range of ±1000 Hz for the first 
satellite is preferably searched. If there is a high degree of 
confidence about the frequency information learned in a first 

40 search, the frequency search range of subsequent satellites 
can be constrained to save time and effort. Alternatively, 
such searching-time savings could be traded for a slower 
clock speed as would be found in less expensive hardware. 
FIG. 10 is a functional block diagram of a second 

45 infrastructure-assisted navigation system embodiment of the 
present invention, and is referred to herein by the general 
reference numeral 1000. The special case in which system 
1000 is needed is where no political or economic opportu- 
nity exists to include a GPS receiver at a cellular-telephone 

50 base site, e.g., system 100 (FIG. 1). So a cell-phone and GPS 
receiver combination permanently parked at a stationary 
location can be adapted to serve instead. A business model 
embodiment of the present invention deploys apparently 
ordinary cellphones on fixed station in each cell site where 

55 business coverage is needed but not otherwise available. 
Such non-availability can result from hostile base cell-site 
operators, adverse or cumbersome laws, loo little time, too 
much expense, etc. 
System 1000 includes a service provider 1002 connected 

60 to a data network 1004, e,g., the Internet, A cellular base 
station 1006 simply handles ordinary cellular telephone calls 
between stations with telephone numbers. A visiting, mobile 
cell user 1008 is able to make and receive telephone calls 
through the cellular base station 1006. A business that 

65 operates the service provider 1002 has previously deployed 
a permanent, parked cell station 1010 that is within the 
service area of the cellular base station 1006. It too is 
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accessible to other parts of the system via telephone calls 
made through the cellular base station 1006. 

The service provider 1002 includes a GPS receiver 1012, 
a network server 1014, and a database 1016 for storing 
temperature, time, and frequency data and models. The 5 
cell-site 1006 comprises a cellular transceiver 1018 and no 
GPS receiver. The cell-phone 1008 also includes its own 
GPS receiver 1020 and cellular transceiver 1022. The 
permanent, parked cell station 1010 has a GPS receiver 1024 
and a cellular-phone-scrvice-band transceiver 1026 that is lo 
used to track the frequency and phase of signals it receives 
from the cell-site 1006. The network server 1014 receives 
and distributes lime information and other data from both 
the cell-site and cell-phone GPS receivers over data network 
1004. A database 1016 of positions and their time and 15 
frequency offsets from GPS time is collected for later 
reference. 

As in system 100 (FIG. 1), three hardware-based mea- 
surements can be used to synchronize the local clocks in two 
independent systems, e.g., (1) time-diff'erencc, (2) frequency 20 
difference, and (3) crystal temperature-versus-frequency 
models. Any or all three can be used to synchronize clocks. 
In a timcDiff observable, if the time difference between time 
events from two different time sources is known, the time of 
an event at one time source can be used to compute the time 25 
of the event at the other time source. In a freqDiff 
observable, if the frequency difference between two clocks 
is known, then finding the frequency of one can be used to 
predict the frequency of the other. In a tempMeas 
observable, a temperature-frequency calibration model is 30 
used to predict the frequency of a local GPS clock by 
measuring the crystal temperature. If one of the time sources 
is a GPS receiver, other clocks can be related to GPS time. 
Such potentially allows devices all around the world to be 
time and frequency synchronized to a common and stable 35 
reference, for example, the GPS system atomic clocks. 

Passive network synchronization embodiments of the 
present invention use the GPS receiver 1024 to measure the 
cell-site 1006 clocks in GPS time. Then the cell-phone GPS 
receiver 1020 clock and the cell-phone 1022 clock can be 40 
synchronized to GPS time by the network server over a 
network infrastructure, e.g., the Internet. A typical GPS 
receiver generates a one-pulse-per-second (1-PPS) lime ref- 
erence that can be compared and related with limeDiff 
circuitry to a 1-PPS clock derived from each of the cell-site 45 
and cell-phone cellular-telephony radio carriers. Other 
observation intervals besides one second would also work 
the same way, but the mathematics are simpler when a one 
second time interval is used. 

Embodiments of the present invention can use three 50 
different communication links to transfer time and frequency 
data. One hoks the GPS reference receivers 1024 with the 
network server 1014. A second carries software-API or 
packet-communication between the cell-phone 1022 and 
GPS receiver 1020 in the cell-phone, either directly in the 5S 
cell-phone or indirectly through the network server. Sharing 
data directly between the GPS receivers helps improve 
system performance. A third connection links the network 
server 1014 to both the cell-phone 1022 and GPS receiver 
1020. 60 

Time-difference-of-arrival measurements (TDOA), e.g., 
from a cell-phone 1010 to a cell-phone 1008, can be used for 
positioning inside a combined cell-phone and GPS receiver 
to improve positioning availability and accuracy. Time- 
stamps sent from the cell-phone 1010 to the cell-phone 1008 65 
can be combined with a priori knowledge of the cell-phone 
1010 position, the cell-phone 1010 clock bias and the 
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cell-site 1006 drift from GPS time. Such time-stamps can be 
used to transfer time information from the cell-phone 1010 
to the cell-phone 1008 with an accuracy affected only by the 
radio-signal propagation distance between. 

The corrective frequency offset of the cell-phone carrier- 
frequency synthesizer loop is again used here in this 
embodiment to estimate the GPS receiver frequency ofket. 
Such estimate is used to reduce the frequency uncertainty 
and thus reduce the time-to-first-fix, improve receiver 
sensitivity, and/or reduce the size and power-consumption of 
the hardware. Digital numeric-controlled oscillator (NCO) 
and analog voltage-controlled oscillator (VCO) are both 
commonly used in such carrier-frequency synthesizer loops. 

Any tcmpMeas measurement can be used to further 
improve the frequency-transfer accuracy passing from a 
cell-phone to a GPS receiver, e.g., by building temperature 
models. The freqDiff observation can be used alone in 
frequency transfers without any timeDiff observation. The 
limeDiff observation can also be used without the freqDiff 
observation for combined time and frequency transfers. The 
time transfer accuracy is preferably improved by reducing 
the timeDiff hardware quantization noise with software - 
based filters. Cell-phone observations of the cell-phone 1010 
frequency and phase can be used in dual-mode positioning 
system alternative embodiments of the present invention. 

The term "passive synchronization" as used herein means 
the time and frequency offsets of the clocks in the cell-phone 
1010 and cell-phone are only observed with respect to a GPS 
receiver clock. Calibrations are made later by sending the 
measured offsets through the communication layers, e.g., 
using transfer control protocol/Internet protocol (TCP/IP) 
packets on the Internet. The real clocks in the cell-phone and 
cell-phone 1010 are not controlled. The signals from these 
clocks are merely observed, and the corrective data obtained 
is not in any way fed back to the source to change the clocks. 
Rather, the data is used mathematically at the network server 
and in the GPS receiver software to assist in the GPS- 
satellite search and then in the position calculation. 

Although the present invention has been described in 
terms of the presently preferred embodiments, it is to be 
understood that the disclosure is not to be interpreted as 
limiting. Various alterations and modifications will no doubt 
become apparent to those skilled in the art after having read 
the above disclosure. Accordingly, it is intended that the 
appended claims be interpreted as covering all alterations 
and modifications as fall within the "true" spirit and scope 
of the invention. 

What is claimed is: 

1. A satellite-navigation system, comprising, 

a first satellite-navigation receiver that provides for a first 
position solution; 

a first wireless communications device collocated with the 
first satellite -navigation receiver and connected to 
share at least one of time and frequency measurements 
with the first satellite-navigation receiver; 

a second satelhte-navigatioo receiver for providing an 
accurate time and frequency estimate; 

a second wireless communications device collocated with 
the second satellite-navigation receiver and able to 
make radio contact with the first wireless communica- 
tions device; and 

a webserver for post-processing of said first position 
solution with said accurate time and frequency estimate 
and for extracting a more accurate second position 
solution for the first satellite -navigation receiver, and 
coimected by a data network to the second sateUite- 
navigation receiver and second wireless communica- 
tions device. 
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2. The system of claim 1, wherein: the second satellite - 
navigation receiver and second wireless communications 
device are at a fixed location and a local differential correc- 
tion measurement is obtained and communicated to the 
webserver for providing increased accuracy position solu- 5 
tions for the first satellite -navigation receiver. 

3. The system of claim 1, wherein: 
the webserver provides information to the first satellite - 

navigation receiver that reduces at least one of local- 
oscillator frequency uncertainty and time uncertainty, 
and provides for improved satellite search times and 
sensitivity. 

4. The system of claim 3, wherein: 
the webserver computes time-difference information by 

comparing time measurements between the first and 
second satellite-navigation receiver as communicated 
over a radio contact and said data network. 

5. The system of claim 3, wherein: 

the webserver computes frequency-difference information ^ 
by comparing carrier firequency measurements between 
the first and second wireless communication devices, 
and such measurements are communicated over a radio 
contact and said data network. 

6. A satellite -navigation system, comprising, ^ 
a first satellite -navigation receiver that provides for a user 

position solution; 
a cell-phone collocated with the first satellite-navigation 
receiver and cormected to share time and frequency 
measurements with the first satellite-navigation 30 
receiver; 

a second satellite-navigation receiver for providing an 
accurate time and frequency estimate; 

a base cellular telephone site collocated with the second 
satellite-navigation receiver and able to make radio 
contact with the cell-phone; and 

a webserver for post -processing of said user position 
solution with said accurate time and frequency estimate 
and for extracting a more accurate second position 
solution for the first satellite-navigation receiver, and ^ 
connected by a data network to the second satellite- 
navigation receiver and base cellular telephone site. 

7. The system of claim 6, wherein: 
the second satellite -navigation receiver and base cellular 

telephone site are at a fixed location and a local 
differential correction measurement is obtained and 
communicated to the webserver for providing 
increased accuracy position solutions for the first 
satellite-navigation receiver. 

8. The system of claim 6, wherein: 
the webserver provides information to the first sateUite- 

navigation receiver that reduces at least one of local- 
oscUlator frequency uncertainty and time uncertainty, 
and provides for improved satellite search times and 
sensitivity. 

9. The system of claim 8, wherein: 
the webserver computes time-difference information by 

comparing time measurements between the first and 
second satellite-navigation receiver as communicated go 
over a cellular radio contact and said data network. 

10. The system of claim 8, wherein: 
the webserver computes frequency-difference information 

by comparing carrier frequency measurements between 
cellular telephone devices, and such measurements are 65 
communicated over a cellular radio contact and said 
data network. 
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11. The system of claim 3, wherein: 

the webserver collects temperature measurement 
(tempMeas) observables and stores them in a database 
from which it later builds a temperature -frequency 
model tised to improve receiver sensitivity and acqui- 
sition performance in a software-compensated-crystal- 
oscillator (SCXO). 

12. The system of claim 1, wherein: 

a minimum of one GPS satellite and three cell. site base 
stations are used in a dual mode system configuration. 

13. The system of claim 3, wherein: 

the webserver includes a global positioning system (GPS) 
digital signal processor (DSP) application specific inte- 
grated circuit (ASIC) that integrates time-difference, 
frequency-difference, and temperature -measurement 
infrastructure-information for improved satellite acqui- 
sition and receiver sensitivity. 

14. The system of claim 1, wherein: 

position-solution calculations are distributed amongst the 
webserver and the wireless communications devices. 

15. A satellite-navigation system, comprising, 

a first satellite-navigation receiver that provides for a first 
position solution; 

a first cell-phone collocated with the first satellite- 
navigation receiver and connected to make and share at 
least one of time and frequency measurements made 
with the aid of the first satellite-navigation receiver; 

a second satellite-navigation receiver for providing an 
accurate time and frequency estimate; 

a second cell-phone collocated with the second sateUite- 
navigation receiver and able to make a telephone 
connection with the first cell-phone through a cellular 
base station located in the area; and 

a webserver with telephone access to each of the cell- 
phones via the cellular base station, and for post- 
processing of said first position solution with said 
accurate time and frequency estimate and for extracting 
a more accurate second position solution for the first 
satellite- navigation receiver, and connected by a data 
network to the second satellite-navigation receiver and 
second cell-phone; 

wherein, the first satellite-navigation receiver and first 
cell-phone are permanently parked at a stationary loca- 
tion and are continually accessible to the webserver, 

16. A navigation system comprises a cell-phone, a base 
cellular telephone site, and a webserver, wherein each are 
paired with a global positioning system (GPS) receiver; 

wherein, a GPS receiver associated with the cell-phone is 
aided with information received from the cell-site and 
webserver that help reduce satellite search uncertainty 
with time-difference and/or frequency-difference mea- 
surements taken with data collected using local cell- 
phone and GPS receiver time and frequency; and 

wherein, correction information is used in post -processing 
of any velocity solutions computed by the cell-phone. 

♦ * * * * 
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